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ABSTRACT 

A  method  for  analyzing  the  performance  of  a  ramjet  engine  at  subsonic 
flight  speeds  is  presented.  The  absence  of  a  known  choked  point  (M*l) 
in  the  engine  necessitates  an  iterative  solution.  A  modified  ideal  gas 
analysis  is  utilized.  Considered  are  the  conventional  ramjet  with  liquid 
fuel  injection  and  an  ejector  ramjet  using  vaporized  fuel  injected  into 
the  engine  at  supersonic  velocities.  In  the  latter  case,  the  fuel's 
momentum  is  significant  and  the  ejector  action  draws  additional  air  mass 
into  the  engine,  which  must  be  considered  in  the  analysis.  The  method 
presented  compares  the  two  engine  cycles  at  several  subsonic  flight  speeds 
for  both  JP-4  and  propane  fuel.  The  effects  of  several  levels  of  component 
efficiencies  are  considered. 
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SYMBOLS 


SYMBOL 


EXPLANATION 

area,  ft2 

inlet  cowl  area,  ft2 

burner  drag  coefficient 

thrust  coefficient,  F/q0A3 

thrust  coefficient  with  additive  drag,  F%  A- 

pressure  coefficient  for  diffuser 

net  internal  thrust,  Ibf 

net  internal  thrust  minus  additive  drag,  lbf 

fuel-to-air  mass  ratio 

conversion  factor  =  32.17405  lbm  ft/ lbf  sec2 

fuel  specific  impulse,  1bf/ibm/hr. 

fuel  specific  impulse  with  additive  drag,  1bf/lbm/hr. 

mass  rate  of  flow,  Ibm/sec 

mach  number 

molecular  weight, 

Ibm-mole 

dump  pressure  loss  e.vpone..L 
pressure,  lbf/ft2 
total  pressure,  lbf/ft2 
stream  thrust,  lbf 

free  stream  dynamic  pressure,  !b*7ft2 

universal  gas  constant  =  1545  -b-  ; bf _ 

It^-molfe  -  °R 

temperature,  °R 
total  temperature,  °R 
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SYMBOL 

(ATT)  actual 

( aTt)  Ideal 

u 

X 

Y 

Z 

7 

Vc 

P 


SUBSCRIPTS 
0  -  5 
P 

* 

C 

T 


SYMBOLS  (CONTD) 

EXPLANATION 

combustor  actual  total  temp,  rise,  °R 
combustor  ideal  total  temp,  rise,  CR 
velocity,  ^/sec 
X  functic 

Y  function  See  Section  III 

Z  function 

ratio  of  specific  heats 

combustion  efficiency  =  (aTt)  actual^  ATy)  ideal 
density;  1bm/ft3 


engine  stations  (see  Figures  1  and  2} 
ejector  ramjet  primary  chamber 
ejector  ramjet  primary  nozzle  throat 
ejector  ramjet  primary  nozzle  exit 
total  conditions 
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SECTION  I 
INTRODUCTION 


The  purpose  of  this  report  Is  to  analyze  the  performance  of  two 
ramjet  engine  cycles  operating  at  subsonic  flight  conditions.  One  Is 
a  conventional  liquid-fueled  ramjet  and  the  other  is  an  ejector  ramjet 

v 

that  uses  gaseous  propane.  Both  JP-4  and  liquid  propane  conventional 
ramjets  are  considered.  The  ejector  ramjet  introduces  its  fuel  at 
supersonic  velocities  with  a  momentum  high  enough  that  it  might  increase 
the  cycle  pressure  ratio  and  overall  engine  performance  over  that 
obtainable  in  the  conventional  ramjet.  The  magnitude  *•  this  performance 
is  determined  on  an  ideal  cycle  basis. 


In  this  report,  three  types  of  engines  are  analyzed  (see  Table  I): 
(1)  a  propane-fueled  ejector  ramjet;  (2)  a  propane-fueled  ramjet;  and 
(3)  a  JP-4-fueled  ramjet.  Several  efficiency  levels  for  each  type  of 
engine  are  considered. 


First,  each  engine  is  assumed  to  have  no  internal  losses  such  as 
burner  drag,  combustion  efficiency  losses,  diffuser  losses,  or  friction 
losses.  Thes'*  results  establish  the  basic  trends  and  serve  to  determine 
the  maximum  values  possible  for  the  performance  variables.  This  "no 
loss"  case  is  referred  to  as  the  ideal  case. 


Secondly,  component  efficiencies  are  applied  equally  to  all  three 
engines.  Baseline  values  of  these  component  efficiencies  were  considered 
to  be  representative  state-of-the-art  values  for  a  subsonic,  JP-4 
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fueled  ramjet  engine;  values  used  were  a  burner  drag  coefficient  of  4.0 
and  a  combustion  efficiency  of  SOSL  These  results  establish  practical 
performance  estimates  for  the  JP-4  ramjet  and  compare  them  with  those  of 
propane-fueled  engines  having  equal  component  efficiencies. 

Instead  of  using  a  burner  drag  coefficient,  we  applied  anotht 
efficiency  factor  to  the  ejector  ramjet.  Since  the  ejector  ramjet  uses 
the  dump  into  the  combustion  chamber  as  the  flame-holding  device  and 
does  not  have  a  baffle  type  flameholder,  the  burner  drag  of  4.0  was 
replaced  with  an  estimated  dump  loss  correction  based  on  experimental 
results  (Reference  5).  Since  the  ejector  ramjet  has  a  diffuser  section 
ahead  of  the  dump,  experimental  corrections  (References  3  and  4)  were 
applied  to  the  flow  to  account  for  this  diffuser  loss. 

Most  of  the  above  mentioned  parametrics  were  computed  assuming 
variable  inlet  size  and,  therefore,  do  not  consider  additive  drag 
effects.  Several  select  cases  have  been  corrected  for  additive  drag 
to  demonstrate  the  order  of  magnitude  of  the  additive  drag  correction. 
The  same  exit-to-combustor-area  ratio  of  0.55  is  used  for  all  three 
engines. 

A  computer  program  was  written  to  calculate  the  performance  of  these 
engine  cycles.  The  method  is  described  in  Section  III.  Appendixes  I 
and  II  describe  the  computer  program. 
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TABLE  I 

CYCLES  INVESTIGATED 


CYCLE 

PROPANE-FUELED  EJECTOR  RAMJET 
Cycle  1 
Cycle  2 
Cycle  3 


PROPANE-FUELED  RAMJET 
Cycle  1 
Cycle  2 


JP-4- FUELED  RAMJET 
Cycle  1 
Cycle  2 


COMPONENT  LOSSES 


No  Losses 

CDb  =  4  and  r/c  =  90% 

Cp  =  0.51 ,  Nd  =  0.25  and 
Vc  =  90% 

No  Losses 

CDg  =  4  and  iyc  =  90% 

No  Losses 

CD  =  4  and  f}r  -  90% 
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SECTION  II 

DESCRIPTION  OF  CYCLES 


1 .  THE  RAMJET  CYCLE 

The  conventional  ramjet  engine  has  been  described  and  analyzed  many 
times  in  the  literature,  and  this  baric  treatment  will  not  be  repeated 
here.  The  majority  of  these  treatments  consider  the  ramjet  at  supersonic 
flight  speeds  with  a  choked  exit  nozzle,  which  presents  a  convenient 
station  to  begin  analysis.  The  ramjet  operating  at  subsonic  flight 
speeds,  however,  usually  has  no  choked  station  throughout  the  entire 
engine.  Figure  1  presents  a  schematic  of  such  an  engine  and  defines  the 
engine  station  nomenclature.  A  convergent  nozzle  is  usually  employed 
and,  since  the  internal  flow  is  entirely  subsonic,  pressure  changes  at 
any  station  are  felt  throughout  the  engine.  The  mass  flow  entering  the 
engine  will  adjust  itself,  generally,  so  that  the  static  pressure  at 
the  exit  (Station  5)  is  equal  to  the  ambient  pressure.  There  are  a  few 
cases  at  high  subsonic  flight  speeds  where  Station  5  can  be  choked  with 

p5  >  po- 

For  the  purposes  of  this  analysis,  an  ideal  inlet  was  considered 
so  as  to  facilitate  the  presentation  of  data  in  parametric  form.  Also, 
no  capture  area  was  specified.  Therefore,  ideal  one -dimensional  flow 
is  considered  between  Station  0  and  Station  3.  The  results  thus  present 
the  net  jet  thrust  coefficient.  A  known  capture  area  can  be  imposed  and 
the  data  presented  can  be  corrected  for  additive  drag  between  Station  0 
and  Station  1.  This  was  done  for  selected  cases  to  describe  the  method. 
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Component  losses  can  be  considered  if  desired.  Total  pressure 
losses  doe  to  friction  or  flameholders  in  the  combustor  are  defined  by 
a  burner  drag  coefficient,  to  be  defined  in  the  next  section.  A  combustion 
efficiency  based  on  borne1'  ideal  total  temperature  rise  can  be  specified, 
if  desired.  The  nozzle  is  considered  ideal  and  no  losses  are  defined. 

Fuel  is  considered  to  be  injected  at  room  temperature.  The  fuel  mass  is 
considered  in  the  continuity  equation,  but  its  momentum  is  neglected. 

The  analysis  procedure  was  chosen  to  facilitate  rapid  calculation 
and  convergence  on  digital  computer  facilities.  Component  efficiencies, 
engine  geometries,  fuel/air  ratio, ana  flight  conditions  are  assumed. 

An  initial  estimate  of  the  free  stream  area  (or  mass  flow)  is  made  and 
the  properties  at  each  station  throughout  the  engine,  from  Station  0  to 
Station  5,  are  calculated.  The  static  pressure  at  Station  5  is  compared 
with  the  ambient  static  pressure  of  the  known  flight  condition.  These 
must  be  equal  for  a  practical  solution  (except  for  the  one  exception 
where  the  exit  nozzle  is  choked).  If  these  pressures  do  not  match,  the 
estimated  value  of  A0  is  modified  .and  the  calculations  are  repeated.  When 
this  matching  of  static  pressures  has  been  achieved,  the  engine  per¬ 
formance  parameters  (thrust  coefficient  and  specific  impulse)  are  calculated. 


2.  THE  EJECTOR  RAMJET  CYCLE 

The  analysis  of  the  ejector  ramjet  ic,  similar  to  that  of  the  ramjet 
except  in  the  treatment  of  tho  fuel  addition.  A  schematic  of  the 
ejector  ramjet  is  shown  in  Figure  2.  A  fuel  injector  in  the  form  of  a 
primary  rocket  with  a  C/D  nozzle  is  located  at  Station  1.  Fuel  is 
heated  to  vaporization  and  is  then  introduced  into  the  primary.  This 
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Figure  2.  Ejector  Ramjet  Engine 
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fuel  is  gaseous  at  high  pressure  and  elevated  temperature.  The  fuel 
expands  in  the  C/D  nozzle  at  supersonic  velocities,  with  a  momentum 
significantly  high  to  be  considered  in  the  cycle.  The  fuel  and  air  mix 
between  Stations  1  and  2  without  burning.  The  mixed  stream  then  enters 
the  combustor  through  a  diffuser  and  sudden  dump.  This  area  change  serves 
as  a  flameholding  device  to  sustain  combustion.  The  rocket  primary 
acts  as  an  ejector  and  draws  additional  air  into  the  engine.  The  momentum 
of  the  primary  increases  the  cycle  total  pressure  ratio  over  that  of  a 
conventional  ramjet. 

Component  losses  through  the  engine  can  be  considered  in  two  ways. 

A  burner  drag  coefficient  can  be  applied  to  correct  for  all  pressure 
losses,  as  was  done  for  the  conventional  ramjet.  The  losses  associated 
with  the  diffusion  and  sudden  dump  can  be  considered  separately  as  a 
function  of  the  geometry  with  correlations  to  be  described  in  the  next 
section.  A  combustion  efficiency  based  on  the  burner  ideal  total 
temperature  rise  can  be  specified  if  desired.  Incomplete  mixing  of  the 
fuel  and  the  air  in  the  mixing  tube  will  result  in  less  ejector  action 
and  less  mass  flow  through  the  engine  than  the  ideal  case. 

A  provision  is  incorporated  to  account  for  these  estimated  losses. 

The  method  chosen  is  to  arbitrarily  reduce  the  primary  momentum  by  some 
percentage  to  obtain  a  reduced  pressure  level  at  the  end  of  the  mixing 
tube  (Station  2).  This  method  was  chosen  for  ease  in  computer  programming 
and  results  in  a  converged  air  mass  flow  less  than  that  obtained  from  the 
ideal  ejector  ramjet  but  still  more  than  that  possible  from  the  conventional 
ejector  ramjet  at  the  same  condition.  The  component  losses  can  be  applied  to 
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the  degree  necessary  to  reduce  Ideal  ejector  performance  to  the  level 
corresponding  to  actual  test  data.  Mixing  losses  were  not  considered 
in  the  study  results  presented  in  this  report. 
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SECTION  IJl 
ANALYSIS  PROCEDURES 


1 .  GENERAL 

The  ramjet  and  ejector  ramjet  cycle?  can  be  analyzed  using  constant 
gamma;  ideal  gas  equations.  It  is  a  requirement  for  subsonic  engines  that 
the  pressure  at  the  exit  of  the  nozzle  be  equal  to  the  free  stream  pressure. 
It  is  therefore  the  purpose  of  the  cycle  analysis  to  find  the  air  flow 
through  the  engine  which  will  allow  the  nozzle  exit  pressure  to  match 
the  free  stream  pressure  for  the  given  engine  parameters.  To  calculate 
the  nozzle  exit  pressure,  it  is  necessary  to  calculate  the  Mach  number 
and  pressure  at  each  station  of  the  engine  beginning  at  the  inlet. 

There  is  one  exception  to  the  above  criteria.  If  the  nozzle  is  choked 
the  exit  pressure  can  be  greater  than  or  equal  to  ambient  pressure 

The  Mach  functions  X,  Y,  and  Z  will  be  used  to  simplify  the  analysis 
p-^cedure  in  this  report.  The  following  is  a  brief  discussion  of  these 
functions.  The  equations  of  continuity,  and  momentum  and  energy  are 
used  to  relate  the  conditions  at  one  engine  station  to  another.  The 
continuity  equation  in  its  simplest  form  states  that  the  mass  flow  rate 
at  Station  i  is  equal  to  the  mass  flow  rate  at  Station  i+1.  This  is 
wri tten 

r"i  =  ^i+i  ^  ^ 

Using  the  ideal  gas  law  p=  PM# 

"rt  " 
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and  the  following  relationships 

m  =p  a  V 
Tt_=  I  +  M  2 

_r 

y- 1 


bquation  1  can  be  written 


PT  AM 

VTr 


r- 1 
2 


y+i 
m2)  — 


R 


+  m  added 


This  equation  is  simplified  by  defining 


This  gives  the  final  Torm  of  the  continuity  equation 


|  Pip  AX  /  9c  M  vf 

* 

ft-  A X  f~\&* 

+  sii  odd  3d  (2) 

Ktt  J  * 

1 

f 

</Tt  V  R 

•  *< 

i-f  1 

The  momentum  equation  for  a  fricticnless,  consvar.fc  duct  is 


+  PA 
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equation  oar,  be  expander!  In  a  manner  similar  to  the  continuity 


equation  to  give 


[?T  AZ]i  s[PTAZ]<+I 


where  l  is  defined  as 


.+y«2  j: 

**oTEuf»r-i 


The  quantity  PTA2  1«  called  stream  thrust  and  Is  denoted  by  the  symbol  P. 
If  the. continuity  and  momentum  equations  are  combined,  the  following 

equation  is  obtained 

_  (4) 

ss  AZ!L  JJL-Z. 


TTMv 


where 


Equations  2.  3.  and  4  are  used  to  calculate  the  flow  conditions  through 
the  inlet.  A  complete  explanation  of  these  equations  is  given  In 
Reference  2.  The  Mach  number  functions  X,  V.  and  I  are  used  to  determ, oe 

the  Mach  number. 

in  this  study,  all  engine  areas  and  the  fuel/air  ratio  are  assumed 
hoowo,  as  war  as  the  free-strcam  temperature,  pressure,  and  Mach  number. 
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The  Ideal  temperature  rise  in  the  combustor,  the  molecular  weight,  and 
the  gamma  at  the  exit  of  the  combustor  are  tabulated  as  a  function  *-  the 
air  total  temperature  and  the  fuel/air  ratio.  Values  of  combust' 
efficiency,  .,c»  and  burner  drag  coefficient,  CDB,  were  assumed  to  pent-.it 
us  to  evaluate  their  effect  on  engine  performance. 


The  combustion  efficiency  is  defined  as 


7)c  - 


ideoi 


The  burner  drag  coefficient  is  defined  as 


Ps  -  P4 


08  ~  in  a3  m*  p3 

where  the  term  P3  -  P4  is  the  drag  loss  caused  by  the  flarnaholder. 


2.  RAMJET  ANALYSIS 

The  ramjet  analysis  consisted  of  the  following:  First  we  assumed  an 
inlet  capture  area.  Then  we  used  the  continuity  equation  to  calculate  the 
conditions  at  the  entrance  to  the  combustor,  and  the  momentum  and 
continuity  equations  to  calculate  the  conditions  at  the  exit  of  the 
combustor.  Finally,  we  used  the  continuity  equation  to  relate  the 
condition  at  the  exit  of  the  combustor  to  the  conditions  at  the  exit 
of  the  nozzle. 


»ao»  s 
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The  method  of  analysis  is  as  follows: 

(a)  Assume  an  A0. 

(b)  From  continuity,  calculate  the  Mach  number  at  Station  3. 
Assuming  the  inlet  isentropi cally  diffuses  the  air  and  the  value  of 
gamma  at  both  stations  is  1.4,  then 


X3 

Mg 


r  X, 


A, 


=  MX,) 


(c)  Calculate  the  free-stream  total  pressure. 

Xo 


PTO  =  po  M  “T-1  ^  )y°J 


Since  isentropic  diffusion  has  been  assumed  in  the  inlet,  it  follows 
that 

PT3=  pT0 

(d)  Calculate  the  mass  flow  rate  of  air 

=  Po  A u  Mo Mwo 
"  R  Ty 

(e)  From  the  given  fuel /ai r  ratio,  calculate  the  fuel  flow  rate. 


P*f  =  mA  f/o 
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Once  the  pressures  have  be  .i  satisfactorily  matched,  the  engine  performance 
can  be  calculated. 

(k)  Thrust  =  PT4  A#  Z #  —  pTo  A0  Z°  ~  po  (  AB~  Ao) 

Isp  =  thrust/rtif 

Note:  If  the  Mach  number  at  any  station  exceeds  one,  reduce  A0  and 
return  to  Step  (b).  If  the  Mach  number  at  Station  5  equals  one  and 
P5  >  P0,  this  is  a  solution. 

3.  EJECTOR  RAMJET  ANALYSIS 

The  ejector  ramjet  was  analysed  in  a  manner  similar  to  the  ramjet. 

In  addition  to  burner  drag  and  combustion  efficiency,  ejector 
effectiveness,  diffuser  and  dump  loss  were  also  considered.  Integration 
of  the  ejector  into  the  ramjet  cycle  analysis  is  the  only  major  deviation 
from  the  previous  analysis. 

The  method  of  analyzing  the  ejector  ramjet  is  as  fellows: 

(a)  Assume  a  value  for  A0. 

(b)  From  the  continuity  equation,  calculate  the  Mach  number  at 
Station  1,  assuming  the  inlet  isentropically  diffuses  the  air,  and  the 
value  of  gamma  at  both  stations  is  1.4. 

X,  =  X0 

A! 

M,  =  f  (X,) 

(c)  Calculate  the  rree  stream  total  pressure 

Yo 

pTo  =Po  it*  ^-r~  )  r°~' 
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Since  Isentropic  diffusion  has  been  assumed  in  the  inlet,  it  follows 

‘hat  PT|  .  pTo 

(d)  Calculate  the  air  stream  thrust  at  Station  1  fiom  the  momentum 
equation. 

Pa,  =  PTo  A,  Z, 

(e)  Calculate  the  mass  flow  rate  of  air. 

=  Po  Ao  Mo  9e  M  w~ 

R  To 

(f)  From  the  given  fuel/air  ratio,  calculate  the  fuel  flow  rate. 

mf  =  mA  t/o 

(g)  Calculate  the  total  pressure  of  the  fuel  injector  assuming 
a  choked  throat,  in  this  analysis  the  ejector  geometry  is  fixed; 
therefore, the  ejector  total  pressure  is  varied  to  match  the  flow  rate. 


where  A*  =  area  of  injector  throat 
X*  =  X  function  at  throat 
Iff  =  total  temperature  of  fuel 


(h)  Calculate  the  Mach  number  at  the  exit  of  the  injector. 
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(i)  Calculate  the  stream  thrust  at  the  ejector  exit., 

=  PTp  Ae  Ze 

(j)  Calculate  the  temperature,  molecular  weight,  and  gamma  of  the 
mixed  fuel  and  air  at  Station  2  by  mass  averaging  the  individual 
properties.  The  individual  physical  properties  are  obtained  from  tables 
of  data  or  empirical  correlations. 

(<  Calculate  the  stream  thrust  at  Station  2  from  the  momentum 
equation.  For  the  Ideal  case 

p£  s  B  +  P« 

If  it  is  desired  to  account  for  inefficiency  of  the  ejector  action,  a 
component  efficiency  can  be  incorporated.  For  example 

\  s  P.  +  V  pe 

(1)  Calculate  the  Mach  number  at  Station  2,  from  the  continuity 
and  momentum  equations. 


Qc  M* 
"2 

m2  =  f  (Y2  ) 
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(n)  Calculate  the  static  pressure  at  Station  2. 

Zi- 

p2  *  pT.(i+  mi  i  '-n 
2 

(o)  Calculate  the  static  pressure  at  Station  2'  if  diffuser  losses 
are  to  be  considered 

P2  =  H  c%r,  m|  + 1  )  p2 

where  Cp  is  an  experimentally  determined  diffuser  loss  factor  obtained 
from  References  3  and  4.  If  icentropic  diffusion  is  assumed 

K  -  p2 

(p)  Calculate  the  Mach  number  at  Station  2'  using  the  continuity 
equation. 


(q)  Calculate  the  total  pressure  after  the  dump  at  Station  3,  by 
using  the  following  equation: 


where  Nq  is  obtained  from  Reference  5. 
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(r)  Calculate  the  Mach  number  after  the  dump  at  Station  3. 


XS 


.  C*k+W"  TT2 
Pj5  A3  9c  M** 


Ms  =  Mxs) 

(s)  Calculate  the  static  pressure  at  Station  3. 

y2 


P3  =  PT3  (1+ 


r2-i 


Mf  )  I*  * 


(t)  Calculate  the  Mach  number  at  Station  4  from  the  continuity 
and  momentum  equat  jns. 


Y4  =  Y3  y^fH.  lc(TU-  TT2)j 
Mw4  L  TT2 


i+  r  •«: 


Tz  [l+  y  Mf(!-Y2  cDB )] 


M4  =  f  {  y4) 


(u)  Calculate  the  total  pressure  at  Station  4  from  the  momentum 
equation. 


rT4 


PT»  [It  -+  C°» 


(v)  Calculate  the  Mach  number  at  Station  5  from  the  continuity 
equation.  It  is  assumed  that  the  values  of  gamma,  total  temperature, 
and  total  pressure  at  Station  5  are  the  same  as  those  at  Station  4. 


M5  =  f  (X5  ) 
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(w)  Calculate  the  static  pressure  at  Station  5. 

r  i  r±- 

(x)  Compare  the  static  pressure  at  Station  5  with  the  ambient 

pressure  P0.  If  these  pressures  do  not  compare  reasonably  well,  adjust 
A0  and  return  to  Step  (b). 


(y)  If  the  static  pressure  and  ?Q  match. 
Thrust  =  pat 


Isp  =  thrust/rru 


T4  A*  Z*  ~  PTo  A°  2°" 


Pq  An  ) 


Note.-  If  the  Mach  number  at  any  station  exceeds  one,  reduce  A„  and 
return  to  Step  (b).  If  the  Hath  number  at  Station  5  equals  one  and 
p5  i  po>  this  is  a  solution. 
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SECTION  IV 
STUDY  RESULTS 

1 .  IDEAL  PROPANE  EJECTOR  RAMJET 

Figure  3  presents  the  parametric  performance  dat.-?  for  a  propane 
fueled  subsonic  ejector  ramjet  at  an  altitude  of  23,000  feet,  A5/A3  =  0.55, 
and  100%  efficiencies.  Plotted  is  the  thrust  coefficients  (Cp)  based  on 
free  stream  dynamic  pressure  and  combustor  area  versus  fuel  specific 
impulse  (ISP).  The  dashed  lines  represent  constant  values  of  fuel-to- 
air  ratio  and  the  solid  lines  represent  constant  values  of  free  stream 
Mach  number.  Several  factors  are  evident  from  this  figure.  First, 
it  is  noted  that  for  this  ideal  case,  as  the  fuel -air  ratio  decreases 
the  fuel  specific  impulse  continues  to  increase  while  the  thrust  decreases. 
Obviously,  the  specific  impulse  must  maximize  at  some  f/a  ratio  and  then 
decrease  as  f/a  ratio  is  lowered  further.  This  will  be  evident  when 
component  efficiencies  are  introduced  into  the  cycle.  The  second 
prominent  feature  occurs  above  the  stoichiometric  fuel /air  ratio  (f/a 
0.064).  As  more  fuel  is  added  above  the  stoichiometric  point  the 
thrust  continues  to  increase.  This  phenomenon  is  not  present  in  the 
conventional  ramjet  because  the  contribution  of  fuel  momentum  is  not 
considered  in  the  ramjet  cycle.  In  the  ejector  ramjet  cycle  as  the  fuel 
flow  rate  continues  to  increase  the  fuel  momentum  increases  and  thrust 
benefits  accrue,  at  a  loss  in  specific  impulse.  Also,  it  is  noted 
that  specific  impulse  improves  sign. ncantly  with  Mach  Number  over  the 
range  considered. 
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2.  PROPANE  EJECTOR  RAMJET  (CDB  =  4.0,  7}c  *  0.90) 

Figure  4  repeats  the  results  of  Figure  3  for  a  propane  ejector  ramjet 
except  that  a  burner  drag  coefficient  of  4.0  and  a  combustion  efficiency 
of  0.90  has  been  included.  If  the  figures  are  compared,  it  can  be  seen 
that  the  performance,  both  thrust  coefficient  and  specific  impulse, 
have  been  lowered  by  including  the  efficiencies.  Also  as  the  fuel-air 
ratio  is  decreased  the  specific  impulse  does  not  continue  to  increase 
as  it  did  in  the  ideal  case,  but  maximizes  between  f/a  '=0.02  and  0.03 
and  decreases  as  f/a  approaches  zero.  Figure  5  considers  additive 
drag  for  an  engine  with  Ac/A3  =  0.2047.  The  design  point  at  which  this 
Ac/A3  was  chosen  is  Mq  =  0.95  and  Cp  =  0.5.  Figure  6  is  a  composite  of 
several  constant  Mach  number  lines  taken  from  Figures  4  and  5.  The 
dashed  lines  in  Figure  6  stop  a.  the  line  representing  full  inlet 
capture.  The  additive  drag  effects  on  engine  performance  are  small  in 
magnitude  but  increase  with  increasing  Mach  number. 

3.  PROPANE  EJECTOR  RAMJET  (DIFFUSER  AND  DUMP  LOSSES) 

It  was  pointed  out  previously  that  the  ejector  ramjet  had  a  sudden 
dump  into  the  combustor,  which  served  as  a  flameholding  device; 
therefore,  perhaps  the  burner  drag  coefficient  of  4.0  used  previously 
was  not  appropriate.  So  in  an  attempt  to  use  component  efficiencies 
consistent  with  the  ejector  ramjet  geometry,  experimental  data  was 
obtained  to  account  for  the  dump  loss  into  the  combustor  and  other  data 
applied  to  the  diffuser  directly  ahead  of  the  dump.  The  method  of 
accounting  for  these  effects  is  described  in  Section  III.  A  diffuser 
loss  factor  Cp  of  0.51  and  a  dump  loss  factor  Ng  of  0.25  were  used 
instead  of  a  burner  drag  coefficient.  A  combustion  efficiency  of  90% 
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was  maintained.  It  can  be  seen  from  Figure  7  that  this  case  with  dump 
and  diffuser  loss  is  worse  than  the  previous  case  with  combustion 
efficiency  of  90 %  and  a  burner  drag  coefficient  of  4.0.  The  absolute 
values  of  Cp  and  ND  are  for  the  specific  geometry  evaluated  and  apply 
only  to  that  particular  case. 

Figures  8  and  9  show  the  effects  of  additive  drag  with  an  Ac/A3  =  0.1825 
which  was  chosen  at  the  design  point  corresponding  to  M0  =  0.05  and 
Cc  =0.5.  As  was  noted  foi  the  previous  case  with  CDg  =  4,  the  additive 
drag  effects  for  this  case  are  small. 

4.  IDEAL  PROPANE  RAMJET 

Figure  10  presents  the  results  for  a  propane-fueled  ramjet  with  no 
internal  flow  losses  and  an  exit  area  ratio  A5/A3  of  0.55  at  23,000  feet 
altitude.  Figure  10  is  a  plot  of  the  thrust  coefficient  versus  fuel 
specific  Impulse.  The  same  trends  are  present  as  for  the  ejector  ramjet 
(Figure  3)  except  that  above  the  stoichiometric  fuel -air  ratio  there 
is  nc  additional  increase  in  the  thrust  coefficient  since  fuel  momentum 
is  not  considered  in  the  ramjet  cycle.  Thrust  coefficient  increases 
with  f/a  ratio  as  specific  impulse  decreases.  Here,  also,  there  is  no 
maximization  of  the  specific  impulse  as  the  f/a  ratio  decreases.  Again, 
the  specific  impulse  increases  as  the  subsonic  flight  speed  increases. 

In  general,  the  ideal  ejector  ramjet  has  better  performance  at  fuel/ 
air  ratios  above  approximately  0.025;  however,  this  comparison  assumes 
no  internal  flow  losses  for  either  engine. 
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5.  PROPANE  RAMJET  (CDg  =  4.0  and  lc  =  0.90) 

Figure  11  presents  the  results  when  using  a  burner  drag  coefficient 
of  4.0  and  a  combustion  efficiency  of  0.90  for  the  propane  ramjet. 

They  are  directly  comparable  to  the  ideal  case.  Figure  10.  While  the 
values  of  thrust  and  specific  impulse  are  lower  for  the  case  with  the 
burner  drag  and  combustion  efficiency,  the  basic  trends  are  the  same 
with  the  exception  of  the  lower  fuel-air  ratios.  As  the  fuel-air 
ratio  decreases,  the  specific  impulse  does  not  continue  to  increase; 
it  maximizes  and  then  decreases  as  the  f/a  ratio  approaches  zero. 

6.  IDEAL  JP-4  RAMJET 

Figure  12  presents  the  ideal  performance  for  a  JP-4  fueled  engine. 

This  data  is  similar  to  the  data  shown  in  Figures  10  and  11  for  the 
propane  ramjet.  Again  no  internal  losses  are  assumed  and  A5/A3  is 
0.55  at  an  altitude  of  23,000  feet.  The  same  trends  are  evident 
although  the  propane  ramjet  has  a  slight  advantage  at  the  low  and  medium 
fuel-air  ratios.  Above  the  stoichiometric  fuel-air  ratio  (0.068  for 
JP-4  and  0.064  for  propane)  the  performance  is  almost  identical. 

7.  JP-4  RAMJET  (C0B  =  4.0  and  ^  =  0.90) 

Figure  13  presents  the  results  for  a  JP-4  ramjet  with  a  burner 
drag  coefficient  of  4.0  and  a  combustion  efficiency  of  0.90.  Figure  13 
is  directly  comparable  to  Figure  12. 

Figures  14  and  15  show  the  effects  of  additive  drag  with  Ac/A3  =  0.1976, 
which  was  chosen  for  the  design  case  of  M0  =  0.95  and  Cp  =  0.5. 
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Oolgn  Point:  M0  *0.98,  Cp*0.6(  Ac/A3  *  0.1978 


I  S  P 
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Figure  14. 


JP-4  Fueled  Ramjet  Performance  With 
Cqb  =  4,  and  i)c  =  90% 
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SECTION  V 
COMPARISONS 


The  ramjet  and  ejector  ramjet  performance  parametrics  shown  herein 
can  be  used  for  several  comparison  purposes  only.  The  results  are  valid 
for  the  assumptions  made.  In  general,  an  application  requires  that  an 
engine  operate  o.er  a  wide  envelope  with  fixed  geometry,  which  necessitates 
considering  additive  drag.  In  addition,  the  external  drag  of  the  engine 
nacelle  must  be  included,  as  well  as  any  vehicle/engine  interference 
drag.  Without  these  specific  effects,  the  following  general  comparisons 
can  be  reached. 

Figure  16  compares  che  parametric  performance  •  -  '?;«  Ideal  engines 
at  Mach  0.7,  The  lowest  point  on  each  line  is  for  .-/a  =  0.02  and  the 
highest  point  is  for  f/a  =  0.1.  As  the  fuel/air  ratio  increases,  the 
thrust  increases  at  a  sacrifice  in  specific  impulse.  At  the  vary  low 
fuel/air  ratios  the  performance  is  nearly  identical.  As  the  fuel/air 
ratio  increases,  the  advantages  of  the  ejector  ramjet  become  apparent. 

The  "X"  mark  on  each  line  indicates  a  stoichiometric  fuel/air  ratio. 

Thrust  levels  below  this  mar  .  indicate  lean  engine  operation,  and  those 
above  this  mark  indicate  fuel  Hch  operation. 

The  propane  ramjet  and  the  propane  ejector  ramjet  can  be  compared 
in  many  ways.  As  noticed  on  the  ramjet  curves,  the  thrust  maximizes 
at  about  the  stoichiometric  fuel/air  ratios;  richer  mixtures  are  of  no 
advantage  to  the  ramjet.  Comparing  the  ideal  engines  at  this  stoichiometric 
fuel/air  ratio  indicates  that  the  ideal  ejector  ramjet  has  a  thrust 
advantage  of  18£  and  a  specific  impulse  advantage  of  11 %.  The  ejector 
ramjet  can  increase  thrust  at  a  sacrifice  in  impulse  by  operating  fuel 
rich;  this  is  no  advantage  to  the  ramjet. 
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Let  us  compare  the  ejector  ramjet  operating  at  a  f/a  =  0.1  and  the 
ramjet  operating  at  stoichiometric.  For  this  case,  the  ejector  ramjet 
has  a  thrust  advantage  of  37%  but  a  specific  impulse  that  Is  only  76% 
of  that  possible  with  the  ramjet.  Figure  17  shows  the  same  comparisons 
at  Mach  0.95.  The  same  comparisons  can  be  made  from  Figure  18  for  Cqb  =  4 
and  i]c  =  90%.  With  both  engines  operating  stoichiometrically,  the 
ejector  ramjet  has  a  17%  thrust  advantage  and  a  10%  specific  impulse 
advantage,  slightly  lower  than  for  the  ideal  case.  With  the  ejector 
ramjet  operating  at  f/a  =  0.1  and  the  ramjet  at  stoichiometric,  the 
ejector  ramjet  has  a  35%  thrust  advantage  but  again  at  76%  of  the 

* 

!  ramjets'  specific  impulse. 

| 

i 

i  If  the  ejector  ramjet  has  the  drag  predicted  from  References  4,5, 

j 

and  6,  and  the  ramjet  has  a  =  4  and  ??c  =  0.9,  we  obtain  the 
following  results.  With  both  engines  operating  stoichiometrically,  the 
ejector  ramjet  has  an  8%  thrust  advantage  and  a  5%  specific  impulse 
advantage  over  the  ramjet.  With  the  ejector  ramjet  operating  at  f/a  =  0.1, 
its  thrust  advantage  over  the  ramjet  is  21%  but  its  specific  impulse  is 
only  65%  that  of  the  ramjet.  Similar  comparisons  can  be  made  at  Mach  0.95 
and  23,000  feet  from  Figure  19;  ft  must  be  pointed  out,  however, 'that 
this  comparison  is  made  st  a  maximum  thrust  level  and  at  a  very  low 
specific  impulse  level,  which  gives  the  maximum  potential  advantage  to 
the  ejector  ramjet.  For  a  cruise  application  a  lean  fuel/air  ratio  would 
likely  be  chosen  to  maximize  specific  impulse;  at  a  condition  of  say 
f/a  =  0.3,  the  advantage  of  the  ejector  ramjet  is  considerably  rtduced 
or  even  eliminated.  For  instance,  at  f/a  =  0.3,  the  ramjet  would  produce 
16%  more  thrust  at  10%  higher  specific  impulse.  One  parameter  which 
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is  important  to  the  effectiveness  of  the  ejector  ramjet  is  the  ratio 
of  the  primary  to  the  inlet  air  stream  thrust.  As  this  parameter 
increases,  the  ejector  ramjet  becomes  more  effective  in  Its  pumping  action. 
Figure  20  is  a  plot  of  this  stream  thrust  ratio  versus  fuel/air  ratio  for 
various  flight  mach  numbers.  As  can  be  seen,  this  parameter  increases 
with  increasing  fuel/air  ratio;  therefore,  the  pumping  action  of  the 
ejector  ramjet  will  be  greater  at  the  higher  fuel/air  ratio;*.  This 
effectively  increases  the  amount  of  air  flowing  through  the  engine, 
thus  giving  more  thrust  than  is  possible  with  the  conventional  ramjet 
at  the  higher  fuel /air  ratios. 
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Figure  18.  Comparison  of  the  Engines  With  Efficiencies  at  Hq  =  0.7 
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SECTION  VI 
CONCLUSIONS 

The  potential  performance  ci  the  ejector  ramjet  and  the  conventional 
ramjet  have  been  determined.  At  high  fuel/air  ratios,  the  ejector 
ramjet  has  a  thrust  advantage  over  the  conventional  r^jet.  The 
relative  ranking  of  these  two  engine  systems  can  change  drastically, 
however,  depending  on  the  internal  flow  losses  and  combustion  efficiency 
assumed  in  the  analysis.  In  addition,  the  relative  advantage  changes 
greatly  with  the  fuel/air  ratio  considered.  The  assumptions  of  C^g  =  4.0 
and  ,Tjc  -  0.90  for  the  ramjet  are  considered  as  state-of-the-art  values 
for  JP-fueled  ramjets.  The  ejector  ramjet  losses  assumed  from 
References  3,  4,  and  5  are  considered  representative,  although  data 
from  a  real  engine  of  this  type  is  lacking.  Predictions  of  internal 
drag  in  References  3,  4.  and  5  are  based  on  experimental  data.  Comparing 
these  cases  shows  that  the  ejector  ramjet  has  an  advantage  at  the  high 
fuel /air  ratios  and  the  conventional  ramjet  has  an  advantage  at  the 
low  fuel/air  ratios.  The  reason  for  this  difference  is  that  with  large 
fuel/air  ratios  the  ejector  pumping  action  is  greater  and  the  cycle 
pressure  is  increased,  while  at  the  lower  fuel/air  ratios  the  ejector 
pumping  action  is  less.  This  is  directly  related  to  the  momentum  ratio 
of  tne  ejector  to  the  inlet  air  stream  which  increases  as  the  fuel/air 
ratio  increases. 

The  data  contained  in  this  section  is  parametric,  with  no  fixed 
inlet  size.  A  real  engine  with  a  known  capture  area  will  have  an  actual 
thrust  lower  than  that  estimated  herein  when  additive  drag  and  external 
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drag  are  Included.  This  was  Illustrated  In  Section  IV  for  one  particular 
design  point.  This  thrust  decrement  should  affect  each  engine  similarly, 
however,  and  should  not  change  the  relative  ranking  derived  from  this 
comparison. 
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APPENDIX  I 

ENGINE  PERFORMANCE  COMPUTER  PROGRAM 
INPUT  AND  OUTPUT  PROCEDURES 

For  ease  of  operation,  the  data  read  into  the  program  has  been 
divided  into  four  sets:  (1)  the  fuel  data  set,  which  contains  the 
tables  of  gamma,  molecular  weight,  and  temperature  rise  for  the  combustion 
products  as  a  function  of  initial  temperature  and  fuel-air  ratio;  (2)  the 
engine  geometry  and  the  efficiency  parameters,  initial  values  of  which 
are  built  into  the  program;  since  this  set  of  data  is  entered  in  Namelist 
form,  only  those  parameters  having  values  different  from  the  initial 
valuer-  need  be  entered;  (3)  flight  parameters  at  which  the  engine  is 
to  operate,  including  the  Mach  numbers,  altitudes,  and  fuel -air  ratios 
for  which  engine  performance  is  to  be  calculated;  (4)  the  job  title 
and  the  job  code.  The  order  of  the  first  three  data  sets  in  the  data 
deck  is  not  fixed,  but  the  fourth  set  must  appear  last. 

Each  data  set  is  identified  by  a  key  word  which  alerts  the  program 
that  the  following  data  belongs  to  a  particular  data  set.  The  key  words 
corresponding  to  the  above  four  data  se^  arc:  FUEL,  GEOMETRY,  RANGE, 
and  PROBLEM.  Each  keyword  must  start  in  column  one.  Tables  I  through 
IV  display  the  form  of  all  the  input  data  cards. 

Table  I  shows  the  format  for  the  fuel  card  set.  Card  1  contains 
the  word,  FUEL,  starting  in  column  one.  Nothing  else  appears  on  this 
card.  Card  2  contains  two  numerical  values:  the  number  of  fuel -air 
ratios  to  be  entered  later  in  columns  1-10,  and  the  number  of  initial 
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air  total  temperatures  in  columns  11-20.  Card  3  gives  the  list  of  fuel- 
air  ratios,  starting  in  column  11,  with  six  numbers  per  card;  up  to 
three  cards  may  be  required.  The  first  ten  and  the  last  ten  columns 
of  these  cards  are  reserved  for  identification  data,  (This  identification 
data  is  not  used  by  the  computer.)  The  other  lists  of  data  in  this  set 
are  entered  on  the  same  format.  Each  list  begins  on  a  new  card. 

Table  II  shows  the  variables  that  are  entered  on  a  Namelist  card. 

A  description  of  this  type  of  data  entry  is  given  in  the  Fortran  Extended 
manual . 

Table  III  shows  the  format  for  the  flight  parameters.  Card  1  contains 
the  word,  RANGE,  starting  in  column  one.  The  second  card  contains  the 
number  of  Mach  numbers,  number  of  altitudes,  and  the  number  of  fuel- 
oir  ratios.  Ten  spaces  are  allotted  per  number,  starting  in  column  one. 
The  third  card  contains  the  list  of  Mach  numbers,  where  each  number  is 
allotted  ten  spaces.  The  other  two  lists  are  similar,  except  that  the 
fuel-air  ratio  list  may  require  more  than  one  card  tc  complete  the  list. 
Figure  21  shows  a  typical  data  deck. 

The  printed  output  from  the  program  gives  the  cycle  performance  and 
many  engine  parameters.  Line  1  shows  the  problem  title  and  the  altitude. 
Line  2  shows  the  capture  area  in  square  feet,  the  conventional  thrust 
in  pounds,  the  corresponding  thrust  coefficient,  specific  impulse, 
specific  fuel  consumption,  fuel-air  ratio^and  the  flight  Mach  number. 

Line  4  shows  the  values  of  thrust  in  pounds,  thrust  coefficient,  specific 
impulse,  and  specific  fuel  consumption,  which  have  been  corrected  for 
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additive  drag.  Line  5  presents  the  engine  stations  and  serves  as  a 
title  for  the  data  immediately  below.  Column  titled  E  presents  data  for 
the  exit  of  the  ejector,  which  is  used  only  for  ejector  ramjet  problems. 
Line  6  gives  the  Mach  number  at  each  station.  Line  7  presents  some  of 
the  important  values  of  ganma  that  were  used.  Line  8  shows  the  flow 
area  in  square  feet  at  each  engine  station.  Line  9  shows  the  pressure 
in  atmospheres  at  each  engine  station.  Line  10  shows  the  total  pressure 
in  atmospheres  at  some  of  the  engine  stations.  Line  11  gives  the  total 
temperature  in  °R  at  some  of  the  important  engine  stations.  Line  12 
shows  the  stream  thrust  in  pounds  force  for  some  stations.  Line  13 

shows  the  molecular  weight  at  two  stations.  Finally,  the  last  line  shows  a 
convergence  parameter  titled  cycle,  the  free  stream  pressure  in  lbs/ft2 

the  pressure  at  the  engine  exit  in  lbs/ft2,  the  air  flow  rate  in  lbs/sec, 
and  the  fuel  flow  rate  in  Ibs/sec.  A  sample  output  is  shown  in  Figure  22. 
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TABLE  I  -  FUEL  CARDS 


Card 

Order 

Contents 

Format 

1 

FUF.L 

AID 

2 

Number  of  fuel-air  ratios  (max  value  -  18) 
Number  of  initial  air  temps  (max  value  -  12) 

2110 

3 

List  of  fuel-air  ratios 

10X.6E10.0 

4 

List  of  initial  air  temperatures 

10X,6E10.0 

5 

List  of  temperature  rise  data  corresponding 
to  the  fuel -air  ratios  and  the  initial  air 
temperatures . 

10X,*Ei0.0 

6 

List  of  molecular  weight  data  corresponding 
to  the  fuel-air  ratios  and  the  initial  air 
temperatures . 

10X,6E10.0 

7 

List  of  gammas  corresponding  to  the  fuel-air 
ratios  and  the  initial  air  temperatures. 

10X,6E10.0 
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TABLE  It  -  VARIABLES  IN  GE0M  NAMELIST 

The  ikey  word  GEOMETRY  precedes  the  namelist  data.  This 
word  is  read  in  on  a  A10  format. 


Variable  Type 


A1  R 

ASTAR  R 

AE  R 

A2  R 

A3  R 

A5  R 

DUMPL0S  L 

ETAF2  L 

ETAFE  L 

ETAMIX  R 

TTF  R 

A2P  R 

ND  R 

CPR  R 

Cdb  R 

nc  R 


Value 

Before  Definition  §  Comments 

Read _ 

1.23  Area  of  station  1  in  sq.  ft. 

0.00753  Area  of  ejector  throat,  sq.  ft. 

0.030121  Area  of  the  ejector  exit,  sq.  ft. 

1.2601  Area  of  station  2  in  sq.  ft. 

5.2414  Area  of  station  3  is  sq.  ft. 

2.8852  Area  of  station  5  in  sq.  ft. 

FALSE  Calculate  diffuser  §  dump  losses  if  true 

FALSE  Use  a  fraction,  ETAMIX,  of  the  ideal 
momentum  at  station  2  if  true 

FALSE  Use  a  fraction,  ETAMIX,  of  the  ideal 

ejector  momentum  if  true 

0.0  Mixing  efficiency 

1300.0  Total  temper.'  .  e  of  ejector  flow  in  °R 
2.52  Area  of  station  2’  in  sq.  ft. 

0.25  Dump  loss  parameter 

0.51  Diffuser  performance  parameter 

0.0  Burner  drag  coefficient 

1.0  Combustion  efficienty 
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Card 

Order 

TABLE  3E  -  FLIGHT  PARAMETERS 

Contents 

Format 

1 

RANGE 

AID 

2 

Number  of  Mach  numbers  (Max  -  8) 

Number  of  altitudes  (Max  -  4) 

Number  of  fuel-air  ratios  (Max  -  20) 

3110 

3 

List  of  Mach  number 

8E10.0 

4 

List  of  altitudes 

4E10.0 

5 

List  of  fuel-air  ratios 

TABLES:  -  ENGINE  IDENTIFICATION  DATA 

8E10.0 

Card 

Order 

Contents 

Format 

1 

PR0BLEM 

A10 

2 

Job  title  and  job  code 

(For  an  ejector  ramjet  the  job  code  is 

any  integer  less  than  or  equal  to  0. 

For  a  ramjet  use  any  integer  greater 
than  0.) 

(The  job  title  can  be  any  comment  the  user  wishes 
to  make) 
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Figure  22.  Sample  Output 
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KOMCRJ 


CUC  6600  FVM  V3. Q-251A  QPT  =  1  07/15/71 


PROGRAM  f*Ol  3RJ  (INPUT,  OUTPUT,  T  A“£5=  INR)  f,  T APE6=0UTPUT) 

RrAL  NO 

EXTERNAL  CYCLE,  3UHAC,  AINLET,  Rl NLFT 
LOGICAL  OUNPLOS,  RJ,  TUP,  E  T  A  F  2 ,  P  TAFt 
OIMENSION  AAM?<i),-  AALTU),  APAP(20>,  TITLEU2) 

COMMON  /CYCL/  A5 ,  F5,  RJ,  T TP,  GAM4 ,  WTM4,  T?0,  STAC 
COMMON  /OUMA /  AM«,  COP,  74 P 

COMMON  /A TNL T /  Al,  AO,  AE,  ACT &p.  A  «0 ,  AMI,  AM3CK,  QUMPLOS, 

1  ETA*!);,  ETAF2,  eTAF£,  FAR,  'AMT ,  IH<r,  P!),  PTP,  TC,  TTF 
COMMON/ HUM/  A?,  AT,  AWT,  AM2,FP,r.AM2,Pn,  t  T2,  WA  ,  WF,  A2P,  CPR,  NO 
NAMFLIST/GEO”/  OUMPLOS,  ETft<-'2,  FTAFF,  ETAMIX,  TTF,  Al,  ASTAR, 
t  AE,  A? ,  A?,  A5 ,  A2P,rPP,  NO,  CnR,  FT  AC 
DATA  PUE/4HFUSL/,  GFO/ 8HGEOMETR Y / ,  UNG/5HRAN6S/ , PR0B/7HPP05LEM/ 
ETAF2  =  ,PALSr. 

ETAFE  =  .FALSE. 

TJMPLOS  -  .PfiLSE. 

ttf  =  nco.o 

A?P  =  2.52 
ETAMIX  =  1.0 
C°R  =  C.51 
NO  =  Q,?5 
GOB  =  C.O 
ETAC  =  1.0 

StU1:1...7« 

AP  =  0.03  0121 

A2  -  1,2601  \bei' 

A?  =  6.2414 
Ac  =  2. *352 
AMOLA  =  1,0/28.966 

1  EE  AO  (f ,25  WORD 

2  FORMA'  (A10) 

W-ITE  (6,3)  WORO 

3  FORMAT  (1M  A10) 

IF  (WOPO  .EO.  F'JP)  GO  TO  10 
IF  (WOPO  .60,  GFO)  GO  TO  20 
IF  (WOPO  .EO.  PANG)  GO  TO  111 

IF  (WOPO  .Fq.  PRO!))  GO  TO  701  _ 

WPITE  (6, 25) 

26  FORMAT  (10H  STOP  PROP) 

GO  TO  1000 

10  CALL  FUELOAT (0.0,  C.O,  0-0,  0.0,  0.) 

GO  TO  1 

20  READ  ( 5  ,G  FOH)  ...  .  .... 

GO  TO  1 

111  Rc40  (p,1U0)  NUMAMO,  NUKALT,  NUMrAR 

1110  PORMAT  (3110) 

READ  (6,700)  ( A AMO ( I) ,  1=1, NUMAMO) 

Rp AO  (6,700)  ( A ALT  (  ' ) ,  I  =  1,NUMALT) 

•V  AO  (6,700)  ( A  F  A  R  ( 1 )  ,  _ I  =  1 , NU  M  P  A  ? ) 

GO  TO  1 

700  FORMAT  (8E10.0) 

701  R6-fln  (6,112)  ( TITLE ( I)  »  1  =  1,12),  TCOOF 

112  PORMAY  (12A6,  10) 

IF  ( EOF (5 ) >  1000,  1002 
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MOHCRJ  '  COt;  6600  FTN  V3‘,0-25iA  OPT=L 

"l  U  0  2  COMTIMUr 

RJ  =  .FAL^F,.  _  _ _ 

IP'  tTCOO=  G  T. "  o’i  R  J  =  .TRUE’i . .  '  . 

00  70S  Ji  =  1,  HUHAlT 
ALT  =  AALT(Jl) 

CALL  ATMOS  (ALT,  TO,  0AOSL,  RHO,  TATSL,  P4PSL,  C,  ViS,  l) 
on  =  PAPSL*i44. 0*14.696 
00  70S  J2=1,NUMAM0 
'»mo'=  AA*9(J2>' 
fT 0  =  T Q*TR( 1.4,  AH05 
H° ITt;  (6;  13) 

13  ^ORHAT  <1H0  13( 10H*  *  *  )) 

JT 0  =  P0M1.P  +  0.2*<AMO**2) >**’.S 
IF  (RJ)  GO  TO  IS 

"tthrat  ='  TTF’oTos/i.'o  '  .  . 

CALL  PFOPAN  ( TTHRAT ,  CPF,  HFAKF5 
GAHT  =  r°F/(CPF  -  1.987) 

TC  =  TTF/1.8 

CALL  PPOPAN  (TC,  CPC,  HTr ) 

T*  IR  =  TO/1.8  _ 

'CALL  AIPTKP  (tAiR,“CPA'r HAIR)  '  . . . . 

HT A  r  HAIR  *  (28. 966*(C*AH0>**2)/(2. 0*32.2=778. 15*1. 8) 

IS  00  705  J3  =  l*  HUMFAR 

FAR  =  4FAR(J3> 

CALL  Fl'FLINT  (fTO,  FAR,  OT,  WTM4.  GSM4) 

TT 4  =  TTO  ♦  ET6C*0T 

vF  (RJ)  GO  T0  498  '  *  ‘  . 

AMQLF  =  FAR/4L.0 
TK2  =  350.0 

HT  2  =  (A“OLA*HTA  +  AHOLr*HTF) i (AHOLA  +  AHOLF) 

00  201  1=1,3 

CALL  PPOPAN  ( TK2,  CPFj,  H) 

CALL  4IRTHR~(T><2,  CPA,“H1)' .  . 

CP  =  (AHOLA’CPA  +  AMOLF*CPF) /(AMOLfi  +  AHOLF) 

TK2CX  =  HT2/CP  *  700.0 
TCHEK  =  TK2CK  -  TK2 
201  TK2  =  TK2CK 

TT?  =  TK2*1.8 

~  AMT  =  <1,0  *  FAC1  7  (~ah6l  A  +"  AHOLF) . 

GAH2  =  CP/ (CP  -  1.987) 

498  ACLOW  =  0  •  1* A3 

IF  (.HOT.  RJ)  GO  TO  489 
41  =  A? 

ART  =28,966 

.  TT2  =”  T T 0* . “  '  '  . .  . 

GAH2  =1.4 
489  ACHIGM  =  41 

X5  =  X (GAM4,  1.0) 

X4  =  X5*AS/A3 
AM4  =  XH(X4,  GAH4) 

Z4  =  Z(GA«4,  AH4~)  -  . 

Y4  s  Y  ( GA’14,  AH 4 ) 

Y4P  =  Y4*S0RT(HTH4*TT2/(AHT‘TT4*GAM2) ) 

IF  (RJ)  Y4P  =  Y4P/(1.0*FAR) 

CALL  50LN=W  (BUHAC,  4.0,  AH4 ,  TRAP,  Yl» 
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AM7  =  Yi 

ACT  s  A7*X(1.4, 4N3)/X(1.4t4M9) 

IF  (RJ>  G o  T036 
tn<t  -  o 

ACM  =  A 1 
ACL  =  Al/2.0 

CALL  SOLNEW  (9INLET,  ACH,  ACL,  TRAP,  Yi) 

IF  (.NOT.  TRAJ)  60  TO  460 
WRITE  (6,14) 

FORMAT  (2«HQRINLET  -  TRAR  =  TRUE) 

60  TO  ?'»n 
ACMAX  =  AC 
ACL  =  r.T3*ACT 

IF  (ACL  . 6T .  ACMAX)  ACL  =  0.ts;*ACMAX . . 

INKT  =  6 

CALL  SCLNEKUTNLET,  ACMAX,  ACL,  TRAP,  YD 
IF  (.NOT.  TRAP)  60  TO  17 
WRITE  (6,75) 

FORMAT  (13H  TRAP  IS  TRUE) 

GO  TO  ?41 

IF  (TNKT  .GT.  i)  GC  TO  3? 

IF  (  AR3CK  .GE.  4*3)  GO  TO  37 

ACHIGH  =  ACMAX 

GO  TO  499 

PT3  =  PT0 

27  =  2 (GAM?,  SM7) 

PT4=PT3* ( 23/74-0. 5*COB*GAM2*iH3**2/<2A*RR(SAM2, AMT) ) ) 
f>5  =  PTA/°R(GAM4,1.0) 

IF  ( •  NOT .  RJ )  ACT  =  AC 
IF  (P6  ,Lt.  RO)  60  TO  497 
GALL  CYCLE  (ACT) 

GO  TO  ?4J  .  . 

APHIGH  =  ACT 

CALL  SOLMrW  (CYCLE,  ACLOK,  ACHIGH,  TRAP,  YNEW) 

IF  (.NOT.  TR A°)  GO  TO  240 
WRITE  (6,500) 

FORMAT  (  1H0,  5X,  11HTRAP  =  TRUE) 

20  =  •'(1.4,  AMO.)  _  . .  .  .... 

r0  =  RTO*AC*20 

21  =  7(1. 4, AMI) 

Ft  =  3T0*AJ.*21 

THRUST  =  F5  -  FO  •  P0*(A5-AC) 
thRUS'TH  =  F5  -  Fi  -Pa‘(A5-Al) 

FI  SR  =  THRUST/WF  .  . 

FISPN  =  THRUSTN/-HR 
SFC  =  3600.G/FIS° 

S^CN  =  3600.0/FISRN 

OF  N  =  0,5*1, 4*P0*A3TAMQ**? 

C"  =  THRUST/OEN 

£FH  =  THRUSTN/O^fj  _  _  _  „ .  _ . _  . _ 

WRITE  76,390)  ( T I T L F ( I )  ,  1  =  1,12) ,  ALT,  AG,  THRUST,  CF,  FISP, 
?  SFC,  FAP,  CHO ,  THRUSTN,  CFN,  RTSPH,  SFCN 
FORMAT ( 1HO/1H012A6,15X, 4HALT  =  , FI  0,1/1 H  3X , 2H AO, 12X,5HTHRUST, 

1  11X,  2HCF ,  17X,  3HISP,  12X,  3HSFC,  i4X,  3MFAP.,  12X,  3HAN0/ 

2  1H  3E15.4.7F15.4,F15,2/14H  ADDITIVE  0R2G, 2E15.4,  3r 15.4) 
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705  CONTINUE 

GO  TO  1  " 

10  0  0  STOP  _ _ _  _ 
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FUNCTION  CYCLE  (ROOT) 

LOGICAL  RJ,  NSOL,  OUMPLOS,  "TAFP,  ETAFE 
"COMMON  /CVCL/  A5, >5,  RJ,  TT4,  GAM4,  WTM4,  TTO,  ETAC 
COMMON  /OUMA /  AW3,  COO,  YAP 

COMMON  /AINLT/  Al,  AC,  AE,  ASTAR,  AMO,  AMI,  AH3CX,  OUMPLOS, 

1  ETAMIXt  ETAF2,  ETAFE,  PAR,  GAMT,  INKT ,  PO,  Pro,  TO,  TTF 
COMMON/OUM/  A2, A3, AWT,  A M2 ,F2, G AH? , PT3 , TT2 , H A , WF, A2P, CPR,ND 
1  AC  =  ROOT 

xi  =  a c * x ( i • 4 ,  amo)/ai 

AMI  =  XMl’Xl,  1.4) 

PI  =  PTQ/PR ( 1 .4 , AMI) 

HA  =  PC*AC*AM0*SQRT((1. 4*32.  2*26. 966) /( 1545. 264*TO> ) 

FA1  =  PT0*A1*Z(1.4,  AMI) 

HP  =  HA*F  AP, 

IF  '(  '.NOT.  Rj‘f  GO  TO"  5 
4*2  =  AMI 
p-’  =  FA1 
GO  TO  5? 

5  XSTAR  =  X (GAMT,  1.0) 

PTP  =  (HF*SORT{(TTF*1545.264)/(32.2*44.0)))/{ASTAR*XSTAR) 

PTH  =  FT°/PR(GAMT,  1.0)  . 

IF (o~H.LT .PI)  GO  to  10  i 

XE  =  XSTAR*ASTAR/AE  ! 

AME  =  5MX  (XF,  GAMT) 
pp  =  ptp*AE*Z(GAMT,  AME) 

CO  TO  35 

10  AME  =  HF/(R1*'aF*S0RT(GAmT»32.2*44./(1545.254*T1F)) 
r f  =  P1*AE*(1.0»GAMT*' AMP*«2) ) 

OTP  =  P1*°R ( GAMT ,  AME) 

15  F?  =  FA  1  +  FE 

IF(ETAFP)  F2=  ETAMIX*F2 
IF  (ETAFE)  F2  =  .FA1  ♦  ?TAMIX*FF. 

Y2  =  ( ( WA  +  HP) *S0RT(TT2*i545. 264/(32. 2*AHT) ) )/F2 
50  Am?  =  YM(Y2,  GAM2) 

52  Z?  =  Z(GAM?,AM2> 

X?  =  X  (GAM? ,  AM2) 

IF  (RJ)  GO  TO  56 

IF  (OUMPLOS)  GU  T^  57 

56  X?  =  X2*A2/A-« 

°T  3  =  ( (MA^HF) *SRRT(TT2*1545. 264/(32. 2*AHD )) /CA3»X3> 

IF (RJ)  PT 3= ( HA* SORT (TT2* 1545. 264/ (32. 2* AWT) ) ) /(A3*X3 K 
GO  TO  c8 

57  CALL  DUMP  (Z2,  X3> 

5C  A«3  =  XM  (>T*,  GAM2). 

Y3  =  Y  (C AH2,  AMT) 

Z3  =  7  (GAM? ,  AMT) 

°?  =  PT3/PR(GAM2,AM3) 

°3  A  T  M  =  P3/{ 144 ,0*14.596) 

70  Y4  =Y3*SQRT( AHT*TT4/(HTM4*TT?) ) 

Y4  =  Y4* ( 1.0  *  GAM2*AM3**2) /(l. OF  CAM?* AMT** ?* (1.0-0. 5*C09) ) 

IF  (RJ)  Y  4  =  Y4* (i.O+FAR) 

AM4  =  YM ( Y4,  GAM4) 

X4  =  X ( CAM4,  AM4) 

X5  =  X4*  A  T/A  5 
AM5  =  XM(X5,  GAM4) 
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Tk~=  7  (GAM4,  AM~4)  '•  ~  " . 

_ PT 4  =  PT3* (Z3/Z4  -  0, 5»C0B*GAH2*AH3»* ?./  (Z4»PR (GAM2,  AH3)  >  ) _ 

P5  =  PT4/ ?R( G AM4, AM5) 

Z5  -  7(C-AM4,  _AM5) 

F5  ~  ~  *4>A5*Z5 

CY"  •'  1.0  -  P5/P0  _  __  __  _ 

GO  TO  300  ~  . 

_ ENTRY  HCVClE _ _ _ _ 

PT2A  =_F?>(A2*Z2*2116.224)  . . 

_ P2_f  PT2A/PR (GAM2, AM2)  __  _  _  _  _ 

P4  =pT3/( 21 16. 2 24*  PR( 3 A  M4  »  A  M  *  ) ) 

POATH  =  P0/2116.224 

.  P5ATM  =  P5/2116«"i?4  “  . . 

_ PT  0  A  =  PT  0/2116 .224 _ _ _ 

PI  =  Pl/2116 .224 
PT3A  =  PT3/2116, 224 
DT4*A  ='  PT4/2116.224 
_  IF  (RJ)  GO  TO  200_ 

.  PT®A  =  PTP/2116. 224  '  . .  . . 

_ _ pea  =  PTPA/PR(GAHT,AHE) _ _ 

"200  WRITE  (6,250)  AMO,  AMI,  AM2,  AM3,  AM4,  AM5,  »ME>  GAM2,'GAM4, 

_1_  GANT,  AC,  Ai,  A2,_  A3,  _A5,  _AE,  _P0ATM,  PI,  £2,  °3ATM,  P4,  P5ATM, 

2  p'EA^  pT0A,  PT2  A~»’  -PT3A,~  PT4A,  PTPAj  TTCf,  TT2,~  Tt4»  FAi,  F2,  F5, 

3  FE,  AWT,  WTM4  _ 

250  '  FORMAT  {RHOSTATION,  9X,"  1H0,'  15X'j  IHT,  15X7  ~1H  2,  16X,  1H3,  15X, 

_  __  JL_1H4»  lgx>  1H5 ,  15X,  1HE/  5H  MACH,  F15.4,  6F16.4/  4H  GAM,  3?X,  _ 

. . 2  F16.4,  16X ,  F  i  6  •  4 ,  16X,"'“f‘16^  4/  5H  AREA,  2X,  ‘4E16.'4,  16X, . 

3  2E16.A/  5H  PRES,  2X,  7E16.4/  5H  TOTP,  2X,  F16.4,  15X,  3E16.4, 

4  16X,  F 16 . 4/  5H  T OTT ,  2 X ,  S16.3,  16X,  E16.4,"l6X,  E16.4, 

_5  /  7H  STP£A Mj__16X ,  2F16 . 4,  _32X,  F16.4,  F16.4/  7H  MOL  WT, 

~6*  38X,  F7.3,  28X,  F7. 3) 

_  _ WRITT  (6,252)  CYCLE,  PQ,  P5,  WA,  WF _ _ _ _  _  _ 

'252  FORMAT '{7M"CYCLE=,  ~E13.'6,  5X,  3HP0'-,  E13.6,  5X.~3HP5=',  E13.6,  5X 
1  3H WA= ,  Ell. 4,  5X,  3HWF=,  Ell. 4) 

30  RETURN 
ENO 
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FUNCTION  AINLET (A)" 

_ LOGICAL  OUMPLOS.  ETAF2.  ETAFE. ACLOG _ _ 

COMMON/OUM/ 'A?,  A3,  AMT,  AM2,F2,  GAH2,PT3 ,  IT?,  HA,  WF,  A'2P,  CPRjND 

COMMON  /BUMA/  AM3j(_  C08,  Y4P  _  _ _  ___ 

COMMON  /AINLT/  AI,  AC,  AE,  ASTAR,  AMO,  AMI,  AM3CK,  OUMPLOS, 

1  ETAMIXl  FTAF2,  ETAFE.  FAR.  GANT.  INKT .  PO.  PTO,  TO.  TTF 
ACLOC  =  . FALSE ■ 

_ Gp._I.Q_l _ 

ENTRY  BINLET 

.  ACLOG  =  ,  TRUE.  _  _ _  _ _  _ 

1  AC=ft 

INKT_H_INKT  +  J, _ _  _ _ 

XI  -  AC*X(i.4,AM0) /Ai 

_ AM  1=  XM  ( X  1,1.4' _ _ 

PI  =  PT0/DR(1.4,AM1) 

WA  =  P0*AC*AM0*SQRT<(li.4*32.2*2fi.966)/J1545.264*T0)) 

FA1  =  PT0*Ai*z’(1.4,AMl> 

WF  =  WA’FAR  _ _  _  _  _ _ _ 

XSTAR  =  X ( GA  MT ,  1.0) 

_ _ PTP  =  ( HF*SQRT( (TTF*1545. 264)/ (32.2*44.0) ) ) / ( ASTAR»XSTAR) _ 

pfH  =  PTP/PR(GAMT ,  1.0) 

IF  (PTH._IJ.P1)  GO  TO  3 _ _ _ _  _  _ 

XF  =  XSTAR*ASTAR/AE 

AME  =_SMX(XE,GAMT)  _  ______ 

FF  =  PTP*AE*Z(GAMT,  AME) 

_ QO  TO  4 _ 

3  AMF=WP/ (P1*AE*SQRT (GAMT*32. 2*44. /  (1545 . 264*TTF> > ) 

Fr  P1*AE*(1.0+GAMT*(AME**2)  ) 

4  F2  -  PA1  *  FE 

_  IF(ELAF.2.)  _  F2=  _ETAM_IX*F2 _ _ 

IF  (ETAFE)  F2  =  FA1  ♦  ETAMIX’FE 

_  _  _ Y2  s((WA»HF) *SQRT (TT2* 1545 .254/ ( 32. 2* AWT) ) ) /F2 _ 

IF  (.NO'T.  ACLOG)  GO  TO  10 

YCK_r  Y (G_AM2 ,  ..0) _ __  _  _ _ _ 

AINLET  =  (1.0  -Y2/YCK) 

IF  (([AINLET  .GE.  0.0)  ._ANO,_(INKT  .EQ.  1) )  AINLET  =  0.0  _ 

AINLET  =  AINLET*!. 0F3 

_ _ GO  TO  500 _ _ 

10  AM2  =  YM  ( Y2,  GAM?)  ~  . . 

Z?_  =  .Z(GAJ12,  AM2)  _  _ _ _ 

X2  =  X(GA«2,  AM2) 

IF  (OUMPLOS)  GO  TO  5 _ _ _ _ 

X3  =  X2  *  A  2 / A  3 

PT  3  =  ( (MA+WF)*SQRT(TT2*1545. 264/(32. 2*AHT) ) ) /(A3*X3) _ _ 

GO  TO 

CALL  BUMP _ (Z2j  X3)  _ _  _  _ _  _ _ 

AM3CK  =  XM(X3,  GAM?) 

AINLET  =  (AM3  -  AM3CK)/AM3 

IF  ((AINLET  .GE.  0.0)  .AND.  (INXT  .EO.  1))  AI NLET  =  0.0 

AINLET  =  AINLET*! . 0E9 _ 

500  RETURN 

_  END  _  __  _ _ 
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SUBROUTINE  OUMP  (Z2,  "X3T 

C  COMPUTES  DIFFUSED  AND  DUMP  LOSSES _ 

. . COMMON/OUM/  A2,A3,  AWT,  AM2,  F2  ,'G  AM2  ,PT3,'TT2i  H A ,  WF,  A2P,  CPR,'0H0' 

„  PT2  =  F2/<A?*Z2) _  _  __  _ 

P2  =  °T2/PR(GAM2,AM2) 

P7  =  <G.5»CPR*GAM2*<AM2**2)  f  1.0)*P2  _  _ 

AONE  =  <GAM2-1.0)/2.0 

_ CONE  -  -  <_( (WA+WF)»»2)*1545.?64»TT2) /(  (P3»*2>  *  (A?P»*2)  *  (GAM2»32,2 

i  *AWT)>  ~ 

0..=  S^RTU,0-4V  '(E*CONE) _ _ 

AM3P  =  (-1.0+Q) /(*.  AONE) 

IF( (AM3P.LE«  0.) »  OR« (AM3P.GE.1.) )  AH^P=  (-!♦ 0-Q) /(2t 0»AQNE)  _ 

AM3P  =  SORT { AM3P) 

PT3  ~  P3»QR(GAM?.AM3P)  _ _  _ 

PT3  '=  Pf3*EXP(-0.5*nND»GAH2*AM3P**2) 

X^  =  (WA  +  WF)*SQRT<1545.264*TT?/(32.2*AMT>_)/<PT3*43> 

RETURN 

END 
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FUNCTION  8UMACU)  .  .  . 

_ COMMON  /R"HA/  AM3t  COB.  V4P _ 

COMMON/OU’i/  A2,A3,AWT,AM2,F2,GAM2,PT3,Tr2,WA,HF,A2P,CPR>N0 

. SO  =”A**2 

_ 9U**ACL=  Y>P*XU.0>JiAM2^S.aM.1.0-CO3.«,.0A?./iSJi.?.Ui,Jt(G_AM2-l^a)/2,.0  .  . 

1  *SQ)>  -  A 

_ 8UHAC  «  P,UHAC*1 . QE2 _ 

RFTURN 

END  _  _ _ _  _ 
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SUBROUTINE  A1RTHR  <T,  C»,  HA) 

A  =  6,3.96  .  _  „  .  _ 

9  =  1.762E-3 

C..*.  ~2.636E-7  . .  _ 

CP  =  A  +  B*T  +  C*(T**2) 

.  tp  =  3oo. a _ 

HCON  =  A*TC  ♦  0.5*B*(TC**2)  +  C* CTC**3> >3.0 
HA  =  A*T  *  k,5*J*(r**2)  +  C*(T**3)/3.0  -  H30N 
RETURN 

_  ENO  _ _ _ 


OOO  jo  0.0  0.0 
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SUBROUTINE  ATMOS (Z,TM, SIGMA, RHO , THETA , DELTA .C A, AMU, K) 

CALLING “SEQUENCE  . . 

CALL  ATHOS  "(2  ,  TM ,S1GK A , RHO ,T HST  A,0 EL T A, C A,  AMuTkT 


=  GEOMETRIC  ALTITUDE  <CT> 

=  MOLECULAR  SCALE  TEHPERATURE  (DEGREES  RANKIN) 


ATMOSO 01 
ATMOSO  7  2 
ATMOSO 03 
ATMOSO ' 4 
ATMOSO'S 
ATMOSO 06 
ATMOSO  T7 
ATMOSO  •’8 


RHO  =_J1_NSITY_  LB-SEC^’J-ET**^1*!  OR_SLUGS-FT** 3 
THETA  =  RATIO  OF  TEHPERATURE  TO  THAT  AT  SEA  LEVEL 

OcLTA_=  RATIO  OF  PRESSURE  TO  THAT  AT  SEA  LEVEL  

CA  =  SPEED  OF  S"0UN0  (FT/SEC) 

AMU  =  VISCOSITY  COEFFICIENT  (LB-SEC-FT*»2) _ 

K  =  1  NORMAL,  _ 

=  2  ALTITUDE  GREATER  THAN  300039.  FT., 

=  7  ALTITUDE  NEGATIVE __  _ 


DIMENSION  HPRTMB(11),TMB(11) ,SIGMAB(ll> , ALM ( 1 1) , ARAY ( 11 ,4) 
EQUIVALENCE  (  ARAY  (I,  1»jH  PRIMS  ( 1)  )  ,  (  AR_A  Y  (1,_2  )  ,  TMB  ( 1>  >  , 
f  (ARAY (1,3) , S1GMA9 ( 1) ) , ( ARAY (1, 4) , ALM(1>) 


DATA  ( ( ARAY ( I , J ) ,  J 

X _ 0. _ 

if  360  89. 2 3D  , 

X  .'2020.997  , 

X  1941Q9.480  , 

X  173884.510 

X  259186.350  , 

X _ 295275.590  , 

Y  744488.190  , 

X  524974.  380  , 

X  557742. 78G  , 

X  656167.  80  , 

OATA  0  /  0.0 

X _ S  /  198 

X  AMU?  /  3.7 

X  TMZ  /  518 

K=1 

IF  (Z)  26jL.8j.17  _  _ 

K=3 

GO  TO  13 _ 

IF(Z.GT. 700000.) 
HPRIM=  (  RE  /  (RE  tZ_)  )  * 
DO  10  M=l,ll 
1 F  (  HPRTM-HPR IM3 ( H) 
CONTINUE 

_H=12 _ 

M=M-1 

IF(  ALM(M)  )  14,15,14. 
TM=TM9 ( M) ♦ ALN(M) * ( 
SIGMA=rX°((1.0MQ,/ 
GO  TO  20 


=1,4) ,1=1,11)/ 

_ 518.688  , 

389.988  , 
389.988  , 
508.788  , 

_ 508.288  , 

298.188  , 

_ 298.188  . _ 

406.138  , 

2386.188  , 

2566.188  , 

2836.188 
1874  ,176  /  , 

_.J2 _ L-i 

372998E-07  /  , 

.688  / 


K=X+ 1 
Z. 


1.00000305 
2.9706958E 
3.2665751E 
1. 2117870E 
5. 8677311E 
1. 7329156E 
1.792_8595E 
9.3921519E 
7. 7658593E 
5. 6324877E 
2. 5726771 E 
RE 

_ ?Z  _ 

RHOZ 


00  , 
-01  , 
-02  , 
-03  , 


-0.00356616 

0. 

0.00164592 

0. 


-04  ,  ..-0.00246883 

-05  ,  0. 

-06  , _  0.00219456_ 

-08  ,  0.01097283 

-10_,  0. 00548643 

-10  ,  0. 00274320 

-10  ,  0.001  2024 

/  2. 3 855531E  07 

L  2116.2 _ 

0.0023769 


ATMOS010 
ATMOSO 1 1 
ATMOSO 12 
ATMOS013 
ATMOSO  t  4 
' ATMOS015 
ATMOS016 
ATMOSO 17 
ATMOSO! 8 
A  TMOSO 1 9 
ATMOS. .0 
'  ‘ATMOSO?! 
ATMOSO '2 
ATMOSO  >3 
ATMOS024 
ATMOSe~5 

, _ AXMOSO  6 

,  ATMOSO'V 
,  ATMOSO  2  8 
,  ATMDSn" 9 
,  ATMOSO  0 
,  ATMOSO ‘ 1 
,  ATMOSO  32 
,  ATMOSO  3 
,  ATMOSO  7 4 
,  ATMOSO 
/  ATMOSO  76 
/,  ATMOSO ' 8 
/ 1_  ATMOSO  9 
7  ,  ‘  ATMOS040 
ATMOSO'.l 


) 11,12,10 


HPRIM-HPRIMB(M) 
ALM ( H) ) ) * ( ALOG ( 


TM0(M)/TM) ) ) *SIGMA9 (M) 
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15 


TM=TM(J  (M) 

SIC-MA=S  TGMA8  (M)*EXP  ( ■ 


(Q*(HPPIM-hPRIMO(M))“/TM6(M)  ) 


20 


13 


RHO=RHOZ*SlGMA 
THET A=TM/THZ  _ 

DELTA-SIGMA* THETA . 

CA=49.  02l77*Sf!RT  (_TM)_ 

><MU=AHUZ*SQRT (THETA**3) *  ( (TMZ  +  S)  /  (TM+S)  >' 
RETURN 


ENO 


ATMOS064 

ATHOSOf>5 

'ATH0SQG6 


c 
c 
c 
'  c 
c 

Q 

c 

c 

c_ 

c 

c 

c 

c 

c 

c 


Z=ATKN(X, Y,N,K,XI) 


WHERE, . . 


X  -  TABLE  OF  INDEPENDENT  VARIABLE  VALUES, 
C*1Y  RE  Ab^ENOING  OR  DESCENDING). 

Y  -  TABLE  OF  OEPFNrENT"  VARIABLE  VALUFS.  ' 
N  -  NO.  OF  POINTS  IN  TfGLES  X  AND  Y. 


K  -  OEGRFE  OF  INTERPOLATION  DESIRED. 

XI-  X-VALUF  FOR  WHICH  INTERPC.flTION  IS  DESIRED. 


ATKNG0B6 
ATKNOO  1 7 
ATKN0008 
ATKNOO  99 
ATKN0C3  0 
ATKNOO 11 
ATKN0012 
ATKNOO!  3 
_ATKN001_4 

ATKN01T5 


THF  IK-fFPOLATEO  VaL'JE  IS  RETURNED  AS  THE  FUNCTION  VALUE. 


31  CELLS  OF  ULANK  COMHON  ARE  USED. 


ATKN0016 
ATKNOO’ 7 
ATKNOO! 8 
ATKN001 9 
ATKNOO  0 


C 

C 

. 

DIMENSION  X  ( N ) ,  YlN) 

COMHON  11,  Kl,  LI,  LL,  LU 

COMMON  XX( 13) ,  YY(l3) 

OATA  "KM AX/  12/ 

ATKNOO  1 
ATKNOOi 2 
A  TKNOO"-’  3 
ATKNOO  4 
ATKNOO’5 
ATKN002& 

IF  {  K  .GT.  KMA X  .OR.  K  .LE. 

0  )  GO  TO  3C0  ATKN0027 

c 

ATKNCO  l- 

K1=K*1 

ATKN00’9 

IF  (X(N)-X(l))  100,10,10 

A1KNOO  0 

10 

IP  (XI-X(l))  20,20,30 

ATKNOO  1 

?0 

LL  =  0 

ATKNOO  '2 

GO  TC  200 

A TKNOO '3 

TO 

IE  (X(N)-XI)  40, 40, BO 

A fKNOO  '4 

4  0 

LL=N-K1 

ATKNOO ’B 

GO  TO  200 

ATKNOO  7  6 

BO 

LL- i 

ATKNOO’7 

LU=N 

ATKNOO  8 

60 

IF  (LU-LL-1)  180,180,70 

ATKNOO  9 

TO 

1I=(LL+LU) /2 

ATKNOO  ,0 

IE  (X(LI)-XI)  80,80,90 

6TKN0041 

80 

LL  -u  I 

A  TKNOO '.2 

GO  TO  6  0 

ATKNOO’ 3 

BO 

LU  =  L  I 

ATKN0044 

GO  TO  eo 

A  TKNOO’iB 

100 

IF  (Xl-X(i)t  120,20,20 

A TKNOO 46 

1?0 

IF  (X(N)-XI)  130,40,40 

A  TKNOO  ’*  7 

1  T  0 

LL  =  1 

ATKUn0‘8 

LU  =  N 

A7KN0049 

140 

IP  UU-LL-l)  180,180,150 

ATKNOO  50 

1 B  0 

LT-(LL+LU. /? 

ATKN0051 

IF  ( v<LI) ■ XI)  160,173,170 

ATKNOO  62 

160 

LL-L  I 

ATKNOO  6 3 

GO  .0  1  '<0 

ATKNOO  >4 

170 

LL  -L  T 

A  TKNOOB5 

70 


T3*' 


m  I  njihnmr— 


AFAPL-TR-72-7 


•f  'M 


ion 


ATKN 


cbc  66oo  ft:i  vit.o-25ia  opt=i  07/19/n 


GO  TO  140 

ATKN0096 

lao 

LL=LL- (Kl+t) /Z 

A  r<N00“7 

IF  (LL)  20,^00,190 

ATKN0058 

193 

IF  (LL+K1-N)  200,200,^0 

ATKM0099 

990 

90  210  1=1, K1 

ATKM0090 

I1=LL+I 

ATXNOO  e 1 

XX(I)=X(I1)-XT 

ATKN0092 

210 

YY(I>=Y(tl) 

AKNOOC* 

OO  220  1=1. K 

ATKi-'OO  •"  4 

00  220  J- \ > * 

ATKNOD'9 

5?0 

YY(J  +  l)  =  (l./'(XX<Jfl)-XX(I)  >  >  ♦  f  Y  Yf  1 1  *XX(J+1>  -  Y  Y  ( J+ 1 )  *  XX  ( I)  ) 

ATKNOO'9 

AT<N=YY(<1) 

6TKN009* 

2cr'J3N 

A  fKMCO  •>« 

C 

ATKNuO' 9 

’03 

c  9 1 N  T  1000  ,  K 

ATKMOO  '0 

1 C  3  0 

-09*4T  ( ■'H0K=  ,  1 12 ,  ■,,H  IS  rNC0o,rCT  FOR  -JNCTION  ATKN) 

ATKMOO  ’1 

CALL  SYSTCM( 200.0) 

ATKNQO ’2 

?NO 

ATKNOQ 

71 


£^*25  i 
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JTINE  SOLNEW 


COC  6600  FTN  V3.0-251A  .0PT  =  1  07/15/71 


SUBROUTINE  SOLNEH(FX,  ALOW,  HI,  TRAP,  Y) 

C  SOLVES  THE  FUNCTION  FX  _FOR  THE  VALUE  OF  THE  INOFPENQENT  VARIABLE_X 
C  WHICH  HAKES  THE  VALUE  OF  FX  EQUAL  TO  ZERO,  "f HE  VALUE  OF  'x  HOST  LIE 
C  IN  THE  INTERVAL  BOUNOEO  BY  ALOW  AND  HI.  IF  ANY  FATAL  DIFFICULTY  IS 
C  INOOUNTERFn  IN  THE  SOLUTION  THE  LOGICAL  VARIA3LE  TRAP  IS  SET  EQUAL 
C  TO  T°UE . 

EXTERNAL  PX 

_ LOGICAL  .TRAP,._FI_RST  _ _ _  __  _ 

blMEN"sfON  FUl,  X  ( 4) 

TRAP  =  .PALSF. 

FIRST  =  .TRUE. 

_  TOL  =  _1 . 0E~4 _ _ _ _  _ 

H  =  HI 

_ SHAL  .=  ALOW _ _ _ 

’if  ("low  .Tt.  HI)  GO  TO  17"  '  '  . . 

H  =  ALOW  _ _  _ _ _ 

SHAL  =  HI 

17  X(  1)  =  A!  OW  _  _ 

F  ( 1 )  =  FXIXd’j)  * 

y  =  x(i:  _ 


.LT.  TOL)  GO  TO  83 


28 


32 


IF  (ABS  (F(l) ) 

_  X ( 3)  =  HI.  _  _ 

F(3)  =  FX ( X( 3) ) 

Y  =  X(7) 

IF  (ABS (F(3) / .LT .  TOL)  GO  TO  83 

Z  =  . SI GN ( F  ( 1 F ( 3) )  _ _ 

IF  (  (F ( 1) *7)  .EQ.  0.0)  GO  TO  3? 

WRITE  (6.28) 

PORMAT  ( 6 OH  THE  FUNCTIONS  FOR  THE  END  POINTS  DO  NOT  HAVE  OPPOSITE 
1SIGNS) 


SOL 

1 

SOL 

2 

'sol' 

3 

SOL 

4 

SOL 

5 

SUL 

6 

SOL 

7 

SOL 

9 

SOL 

10 

SOL 

~2 

SOL 

13 

SOL' 

14 

SOL 

15 

SOL 

16 

SOL 

17 

SOL 

18 

SOL 

IB 

SOL 

:>o 

SOL 

’i 

SOL 

22 

SOL 

23 

SOL 

'’!* 

SOL 

25 

SOL 

26 

SOL 

■>7 

TRAP  =  . T  \U£ . 

SOL 

30 

GO  TO  83 

.SOL  U 

X  ( 2 )  =  X(l)  -  F (1) * (X ( 3)  -  X(1))/(F(3)-  F (1 ) ) 

SOL 

{2 

F ( 2)  =  FX:Xt2 )J 

SOL 

T3 

Y  ^  X(2) 

SOL 

34 

IF  (ADS(F (2) )  .LT.  TOL)  GO  TO  83 

SOL 

T5 

DO  69  JK=1 , 9 

SOS 

16 

X1SQ.  =  X(l)**2 

SOL 

37 

X1X2  -  X(l)  -  X ( 2) 

SOL 

38 

XI X 3  =  x(11  -  X (3) 

SOL 

39 

X2X1SD  =  X(2)**2  -  X(l)**2 

SOL 

40 

A  =  X1Y3* (P( 2) -F ( 1) )  -  X1X2* (F (3) -F (1) ) 

SOL 

.1 

A  =  A/(X1X3*X2X1SQ  -  XIX2MX(3>**2  -  X1SQ)> 

SOL 

42 

3  =  ( A*  X?  X1SD  -  F  (  2)  P|1))/X1X2 

_ .SOL _ 

43 

C  =  F (3 )  -  A»(X(3)»*2)  -  C*X<3> 

SOL 

44 

Q  =  SQRT(9**2  4.0*A*C) 

SOL 

4  5 

X ( 4 )  =  (-B+Q) /(2.3*A) 

SOL 

46 

Ir( (X(4) .GT.H) ,OR. (X(4) ,LT .SHAL) )  X (4 > = (-B-Q) /(?. 0 *A ) 

SOL 

4  7 

Y  =  X  ( 4 ) 

SOL 

-.8 

IF  (JK  ,.EQ. .  9)._GO_TC  70 

.  SOL 

4  9 

F(4>  =  FX ( X{ 4) ) 

SOL 

50 

IF  (A0P(P(4))  .LT.  TOL)  GO  TO  83 

SOL 

51 

DO  62  T  =1 , 3 

SOL 

52 

IS  =  1  +  1 

SOL 

•» 

DO  62  J--I1  ,4 

SOL 

"4 

72 


ns 


UTINE  S0LN9H 


ror  6600  FTN  V3.0-251A  OPT=l  07/15/71 


100 

110 

115 

63 

66 

63 

69 

70 


81 

82 

86 

87 


IF  (9(,,  ,t_E.  F  ( J)  )  GO  TO  6° 

F6  =  F { I ) 

XS  =  'X(I) 

9(1)  =  F(J> 

X(I)  =  X  C  J) 
r< J)  =  FS 
X(J)  =  XS 
CONTINUE 

19  (.NOT.  FIRST)'  GO  TO  63 
IF  (  ( F ( 1 )  4  F ( 2 ) )  .LE.  9  < i ) )  GO  TO  100 

XN  =  (?. n»X(2)  *  X  (1) ) /T.O 
FN  =  Fy (XN) 

IF  (  (F ( 1)  4  FN)  .LE,  F(l))  30  TO  110 
F ( ? )  =  FN 

X(?)  =  XN . 

GO  TO  FT 

IF  (<F<7)  4-  F  ( <4  >  )  ,GE.  9(4))  G"  TO  6’ 
XN  =  ( X  ( 4  )  4  2 . 0  *  X  (  3) )  /  3 . 0 
FN  =  FX(XN) 

IF  ((F(4>  4 
9(3)  -  FN 

X  ( 3)  =  XN 
GO  TO  63 

9(1)  s  FfJ 
X(l»  =  XN 
30  TO  63 
F ( 4 )  =  FN 
X ( 4 )  =  XN 

9IRST  '  .FALSE. 

IF  (fiO'  (9(1))  . GT.  ASS (9(4)))  GO  TO  68 

IF  (  (9(3)4 K (4))  .39.  F ( 4 ) )  30  TO  66 

9(1)  =  9(4) 

X(l)  =  X ( 4 ) 

30  TO  6q 

Ic  (  (9  ( 1)  ♦  F{  2)  >  .ur.  9(D)  30  TO  65 

CONTINUE 
Oo  32  1=1,6 

9Y  =  Fy ( Y ) 

IF  (409(FY)  ,  LI  •  TOD  GO  TO  77 
XL  I  I  =  Y  4  0.01*TOL*Y 
FLIT  r  FX(XLIT) 

XO  =  Y  -  XL  I T 
90  =  Fy.  FLIT 
OIF  =  90/XO 

19  (OIF  .NS,  n,3)  r,0  TO  31 

T040  =  .  TRUr . 

G°  TO  »’ 

Y  =  Y  -  c  y/o 19 

C0NTINU9 

WF I T  9  (6,85) 

63  =  3  47  (HH  SOL  NrW  njr  NOT  OONVr,,GF) 

IF40  =  .T=!J9. 

F  F  T  'J  3  N 

rN0 


SOL  55 
SOL  "-6 
SOL  37 
SOL  68 
SOL  6  9 
SOL  60 
SOL  61 
SOL  62 


SOL 

•3 

SOL 

.4 

SOL 

65 

SOL 

>6 

SOL 

67 

SOL 

'  3 

SOL 

■  0 

SOL 

70 

SOL 

’1 

SOL 

’2 

SOL 

’3 

SOL 

*  H 

SOL 

ys 

SOL 

’6 

SOL 

.’7 

SOL 

7  8 

SOL 

79 

SOL 

30 

SOL 

31 

SOL 

37 

SOL 

73 

SOL 

34 

Ill  C'  X 
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ON  X  CPC  6630  FTN  V7.0-251A  OPT=i  07/15/71 


FUNCTION  X!G,  A)  XI 

PMINFS  X  AS  A  FUNCTION  OF  MACH  MUMPER.  3=GAHHA,  A=NACH  NO.  X  2 

l‘=r,"ri.O  .  . . .  "X  1 

=  A*S'V>T<G/(<1.J+0.5»G1*<A**2)  >**((GU.O)  /Gl)>)  X  4 

RETURN  X  5 

FNO  X  6 


7* 


AFAPL-TR-72-7 


:oh  Y 


CDC  6600  FTN  V3.3-251A  OPT=l  07/15/71 


FUNCTION'  Y(G,  A) 

C  DETERMINES  Y  AS  A  FUNCTION  OF  MACH  MUMPER.  G=GAMMA,  A=MACH  NO. 
'  SO  -  A**2 . . 

Y  =  <A/{1.0*G*SQ))*SW(G*(1.9  +  e.5*<G-1.0>*SQ)) 

RETURN 

ENO 


Y  1 

Y  2 

Y  3 

Y  A 

Y  5 

Y  6 


75 
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JTINF  FUcLOAT 


COC  6600  FTN  V3.Q-251A  0PT=1  G7/15/7i 


SUBROUTINE  FUELOAT  (TO,  FAR,  OT,  WTM4 ,  GAM4) 

DINE NSION  TEHS ( 12)  ,  FARS(18?  ,  Z1  <_18 1 1 2 72(16,12),  Z 3J  19,12) 
~READ(5,2>  NEAR,  NT 
FORMAT  (2110) _ 

REAO  (5,5)  (FARS  (I)  ,  I=i  ,*NFAR) 

REAO  (5,5)  (TENS ( I) ,  1=1, NT) 

FnRHA'T  (10X,  6E 10  •  0 ) 

00  _10  _K=liNL  _ _ _ 

'READ' '(5,5”)  (ZlTi,K‘  I  =1,NFAR)  "  ’  . . . 

10  20_  K=1»NT  _ _  _ _ 

READ  (5,5)  ( Z2  ( I, K)  ,  I  =  lj'NFAR) 


00  TO  K= 1 ,NT 

REAO  (5,5)  ( Z3 (I 

,K), 

1=1, NEAR) 

RFTURN 

ENTRY  FUELINT 

OT  =  B'JYIN  (FAR, 

TO, 

FARS, 

NEAP, 

TEMS, 

NT, 

Zl) 

HTN4  =  BUYIN (FAR, 

TO, 

PARS, 

NFAR, 

TEHS, 

NT, 

Z  2) 

GAM4  -•  BUYIN(FAR, 

TO, 

FARS, 

nfar. 

TEMS, 

NT, 

Z3) 

RETURN 

ENO 


j'  vn 
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CON  5'JYIM  .  COC  6500  FTN  V3.3-251A  OPT=l  07/15/7: 

FUNCTION  OUYIN  CXI,  YI,X,'NX,  Y,  NY,  2) 

OIHENSION  XC18)_,  YC12),  2(18,12),  V(18),  U<12) 

00  6  1  =  1, NY . 

00  5  J=1 , NX 

V(J)  =  7 ( J , I ) 

U(I)  =  ATKN  (X,  V,  NX,  1,  XI) 

OUYTN  =  ATKN  (Y,  U,  NY,  1,  Y T ) 

=>ETUPN 
END 


AFAPL-TR-72-7 
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ION  XM 


COC  660  0  FTN  V3.0-251A  OPT=t  '  07/15/71 


"  FUNCTION  XH*( X » 6 ) 

C  DETERMINES  SUBSONIC  HftCH  NUHBER  AS  A  FUNCTION  OF  X.  G=GAMHA 

"foL  =  X*T  .  OF-8  ‘  . ” 

G1  =  G-1.0 
XN  =  0.0 
FI  =  -X 

IF<ABS<F1)  .LE.  TOL)  GO  TO  26 

XM  =  1.0 _ _  _ 

" . f?"  =  SCRT(G/(  (1.0  +  0. 5*Gi)*»{tG*1.0)  /Cl'))  -'  X  . 

Ic  ( ABS  ( FT )  .  LT._  TOL)  GO  T0  26 
SI  =  SIGN ( FI »  M> 

IF  '  (FI  *■  SI)  .EO.  0.0)  GO  TO  13 
PRINT  lOO',  F3 

100 _ FORMAT  (_3_0H  SANE  SIGN  FOR  X  FUNCTIQI  -  F=,  £13.5)  _ 

XM  if  1.0 

GO  TO  26  _ 

13  XM  =-Fl/(FX  -  FI) 

F?  =  XM*SORT(G/((1.04-0.5*Gl*!XH**2>)**<(G*i.0>/Gi>>>  -  X 
A  =  <(F3-Fi)*XH  -  F2  +'  FI)  7  f  XM*  ( 1 . 0-XM)  )' 

B=F3-F1-A 
C  -'Fi' 

RT  =  SORT  (3**2  -  4.0*A*C) 

XM  =  <-B*RTj /(2.¥*Aj 

IF( (XN.GT.1.0)  .OP.  (XM.LT.0.0))  XM= < -8-RT > / ( 2. G* A) 

00  25  1=1,9 

SO  =  XN**  2 

F~  =  ’  X“*SORT(  G/ ("( 1  .*0*0 Vs  *Gi*ST)  *  *  ( (G +1 . 0)  /01 )  TT  •*  X 
IF { ABS I F)  .LT.  TOL)  GO  TO  26 

FP  =  2.0M1.0-SO)*  (SORT  (G/((1.0*3.5*G1*SO)**((G*1.0)/G1)))) 
1  / ( 2 . 0  *  Gl*SO) 

25  XN  =  XM  -  F/fp 

onlNT__50  _  _  _  _ _ 

'  50"  FORMAT  "(2 3H  X"  FUN'0I9~NUf  CONVERGE)'' 

26  RETURN 
END 


XM . 

1 

XN 

2 

XM 

5 

XH 

6 

XM 

7 

XM 

8 

XN 

9 

XM 

’  0 

XM 

l  i 

XM 

12 

XM 

14 

XN 

15 

XM 

1  6 

XN 

17 

XM 

18 

XM 

20 

XN 

21 

XN 

22 

XN 

■*3 

XM 

■*4 

XN 

25 

XM  ?6 
XN  27 
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ION  YH 


CDC  6600  FTN  V3~.0-251A  OP?=i  07/15/71 


FUNCTION  VM(V,G) 

OETSPHINFS  SUBSONIC  MACH  NO.  AS 
TOL=  Y*1.0E-8' 


FUNCTION  OF 


G=GAHMA 


YH  =  0.0 


FI  =  -Y 

IF  fARF  (FI)", LE.  TOL)  G0~T0"27 
YH  =  1.0 

F3  =  (SQRT<5*(1.0<-0,^*6l))>/(1.0  +  G)  -  Y 
IF  (ARR(F3)  .lt.  TOL) __GO  TO  27  _ 

YN  =-?!/( F3  -  FI) . 

SO  =  YH**  2  _ 

F2  =  ( YH/  ( 1.  6+G*SQ)’)  *SQRT  (G*  <1.  0*0.  5*Gi*SQ)  T 
A  -  ( ( F3-F 1>  *YM  -  F2  +  FI) / < YH* ( 1 . 0-YH) ) 


sea 


ION  SMX  '  ■'  roc  6600  FTN  \/3.0-251A  0PT=1  07/15/7' 


FUNCTION  SIX (X*  Gi 

SMX 

1 

C  0£TFPMIsjr$  SU°F.»SONIC  MACH  NO.  AS  A  FUNCTION  OF  X.  G=GAMHA, 

MAX  HACK  =4 

SMX 

DIMENSION  PA(X),  “A(3) 

SMX 

3 

TOL  =  v*1.0F-5 

SMX 

4 

G1  =  G  -  1.0 

s-ix 

5 

X  A  ;  1 )  =  1.0 

SMX 

6 

X  A ( 2 )  =  4.0 

SMX 

7 

00  IX  1=1,3 

SMX 

n 

fa<I)=XA(I)»S0RT(G/<(1,f.5*G1* <XA{I)**2))**<(G*1.)/G1>) ) 

-X 

SMX 

0 

S«X  =  X  A ( I ) 

SMX 

0 

I F ( APS  <  fa ( I)  )  .LT,  TOU  GO'  TO  32 

SMX 

1 1 

IF ( I.CF.2>  XA(3)=XA(1)-FA(1)*(XA(,)-4A(1)>/(FA(2)-FA(1)) 

SMX 

12 

IX 

CONTINUE 

SMX 

l  x 

X1SO  =  XA(1)»»2 

SMX 

•  4 

XIX?  =  XA ( 1 )  -  XA  ( 2) 

SMX 

15 

X1X3  =  X A  ( 1 )  -  XAC) 

SMX 

16 

X?X1S0  =  yA(2)*»2  -  XA(1)**2 

SMX 

17 

A  =  X1X3»(FA(2)  -  F A ( 1 ) )  -  X1X2* (PA  !3) -FA (1) ) 

SMX 

1  8 

a  =  A/<X1X3»X2XIS0  -  XI X2* ( X 4  <  3) *  *2  -  X1S0) > 

SMX 

'  9 

X  =  (A»X?viSD  -  FA  (2)  «•  FA<1))/X1X2 

SMX 

>0 

C  =  F  A  { 3)  -  A»<XA<3)**2>  -  a*Xi(3) 

SMX 

■1 

0  =  SOPT  ( 9**2  -  A .  0*  A*C  ) 

SMX 

’2 

S'-X  =  (-9+0) /<2.Q*A) 

IF(  (5NX.LT. 1.0)  .0®.  (SNX.GT.t,  0))  SNX=(-3-Q)/(2.0*A) 

SMX 

°  3 

00  xt  1=1,5 

SMX 

’5 

. 

SO  =  S»'X**2 

c  =  SMX*SORT<G/(<I. 0*0, 5*G1*S0)**((G*1. 31/31) )>  -  X 

—  •  -  •  — 

SMX 

*■6 

IF  (Aq'-(c,  ,  L  T .  TOL)  00  TO  •»» 

SMX 

•»8 

cc  =  2. 0*  <1.  j-SO)  *  (SORT  <G/<  <1.0*3 .5*G1  *S9>  »♦<  <G*1.  G)  /Gj.)  )  )  )  /  <2. J  ♦ 

SMX 

1 9 

1  3 1  *  SO) 

SMX 

•n 

31 

q-y  =  cmx  -  f/fo 

SMX 

U 

32 

XFT'JSM 

SMX 

'2 

"NO 

SMX 

•  3 
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"  coc  gsqo  ptn  vT.o-astA  opt=i  07/is/;-i 

FUNCTION  OR(G,A> 

C  OETE^HINFS  THE  RATIO  OF  TOTAL  ->RESS'JRE  TO  ST  A  T I  ^  PRSSS'JRF  AT  A 
C  FUNCTION  or  MACH  HO,  ANO  GANNA.  G=GAHMA  .  A=NACH  NO. 

Gi  =  G-l 

=  c.'\  ♦  o.5*Gi*(A**2))**(r,/r,n 
?*'  TORN 
CN  0 


TR 


CDC  6600  FTN  tf 3.)-25lA  OPT=l  07/16/7 


ton 


FUNCTION  TP(G,A) 

C  _nFT6  =  MIN6S  THE  RATIO  OF_  TOTAL  PRESSURE  TO  STATIC  PRESSURE  AS  A 
C*  FUNoTiOH~OF  HACH  'NO'.  AND  GANNA.’  G=GAMMA,  A  =  HACH  NO. 

TR  =  1.0  ♦  n.5*(G-1.0>»(A»*2> 

RETURN 

ENO 


I 
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CORE 


N*p  08.56,31. 


NORMAL 


-l i-L2 — -™£°S. 

FHA  Lim'jcs  nun7A£.  ru>  Tr-.m - —LL.:= _ 


«««*  F«»'TntEnss«r 


-BEsIr3^-*-0^55- 


■ITl— USER—  -»♦ - CM  i  - - 


CYCLE 


000152 


003511 


ikM^LEQ - COMMON — 

6ycl - 


BUM* 

ainlt' 

OUH 

cycl‘~ 

BUMfl 


TTfoioo 

000110 


AINLET" 


DUMP 
SUMAC  ' 

pR0°AN 

airthr 

ATMOS 

ATKN 

’  SOL NEW ~ 
X 
Y 
2 

eueloat 


BUYTN  .  0100 


XM 
YM 
ZM 
SMX 
TR 
PR 

GET PA”  ~ 

sro-t 

SYSTEM? 

TFENuFJ 

INPUTCS 

.JNPUTNS 

KOOERS"” 

KRAKERJ 

OUT°TCI 

ALNLOGE 

EXPE 

SQRTjj 

X>OYE 


004365 


_ 004673 

005005  ' 

005036 

005124 

005171 

_ 0  05  363 

005616 
006265 
006310 
006334 
006360 
10 


aTklT 

OUM 
"OUM 
BUHA 
AINLT  ' 
OUH 
“BUM  A 


000113" 

000135 

oooToo " 
000110 
"000113 
000135 


000135 

,000110 

000113 

000135 


OUM 


000110 

000135 


010147  ' 

_0 10  371 
010570 
_  010766 
011213 

_ 011231 

Oil? 53 

_ 011272 

012646  " 
_  013631 
013710 

_ 014026 

015205 
..  016601 
020325 
_.020421 
020460  ' 

020524 
020546" 
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